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SUMMARY 


This  work  examines  the  simultaneous  use  of  single-gimbal  Variable  Speed  Control 
Moment  Gyroscopes  (VSCMGs)  as  spacecraft  attitude  control  actuators  and  energy 
storage  devices.  The  resulting  theory  is  then  used  as  the  foundation  for  designing  the 
VSCMG  Workbench,  a  simulation  tool  designed  to  meet  the  Air  Force  Research  Lab¬ 
oratory’s  need  for  a  realistic/flexible  computer  simulation  for  conceptual  analysis  and 
hardware-in-the-loop  testing.  This  tool  allows  both  Georgia  Tech  and  AFRL’s  Space 
Vehicles  Directorate  a  low-cost  alternative  for  analyzing  the  feasibility  of  employing 
a  combined  attitude  control/energy  storage  system  as  well  as  the  technical  details  to 
create  such  a  system  for  different  types  of  spacecraft.  Its  modularity  permits  adding 
more  model  fidelity  in  the  future  with  little  user  training  required. 
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Bodies  : 

P  =  spacecraft  platform  without  VSCMGs 
G  =  gimbal  structure 
W  =  wheel 
Mass  Center  points  : 

P*  =  center  of  spacecraft  body  without  VSCMGs 
G*  =  center  of  VSCMG  gimbal  structure 

W*  =  center  of  VSCMG  wheel  (coincides  with  G*  by  assumption) 

O  =  combined  system  center  of  mass  (assumed  constant  in  frame  J\f  described  below) 

Coordinate  Reference  Frames  : 

B  =  ref  frame  that  rotates  with  the  spacecraft  body  and  is  centered  at  P* 

G  =  ref  frame  that  rotates  with  VSCMG  about  gimbal  axis 
J\f  =  inertial  reference  frame 

Go  =  body  fixed,  VSCMG  installation  ref  frame  (note  that  Q  coincides  with  Go  at  the  initial  time) 


9s,9t,9g  —  Unit  vectors  for  VSCMG  spin,  transverse,  and  gimbal  axes 

9s,9t,9g  =  VSCMG  unit  vecs  expressed  as  components  in  B 

rT  =  torque  required  for  stable  regulation  or  tracking 

Lbn  —  Transformation  from  Af  to  B  —  L^B 

Lrn  =  Transformation  from  Af  to  7Z  =  LTNR 

Lbr  —  Transformation  from  K  to  B  =  L^B  —  LBnLrN 

j3  =  column  matrix  of  euler  parameters  representing  rotation  from  J\f  to  B 

(3r  =  column  matrix  of  euler  parameters  euler  parameter  representing  rotation  M  to 

7Z  where  71  is  fixed  to  J\f  in  the  regulation  problem 

oj  =  body  angular  velocity  with  respect  to  M  written  in  B  (in  expanded  notation  this 
is  written  as  (ojb/n)b) 

cor  =  reference  body  angular  velocity  with  respect  to  A f  written  in  B  for  regulation 
and  71  for  tracking 

*=&u 

d)r  —  ~uir  (zero  for  the  regulation  case) 

Lbn  =  Inertial  derivative  of  LBn  =  j^LBN 
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R-Bii  =  Inertial  derivative  of  Lbr  =  ^ LBr 

8/3  —  error  euler  paramter  column  matrix  representing  a  transform  from  B  to  71 
8oj  —  angular  velocity  error 
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CHAPTER  1 

INTRODUCTION 


1 . 1  Overview 

1.1.1  Problem  Background 

Space  vehicle  programs  consistently  seek  to  reduce  satellite  bus  mass  to  increase 
payload  capacity  and/or  reduce  launch  and  fabrication  costs.  In  addition,  satellite 
system  performance  demands  continually  challenge  space  vehicle  designers.  Specifi¬ 
cally,  larger  space  structures  require  finer,  more  accurate  three-axis  attitude  control 
methods.  One  of  the  most  popular  methods  for  this  type  of  control  is  by  employ¬ 
ing  gas  jet  thrusters  to  provide  the  necessary  vehicle  control  torque.  Unfortunately, 
as  mentioned  in  reference  [1],  the  plumes  of  such  thrusters  often  impinge  on  criti¬ 
cal  vehicle  components  such  as  communications  equipment,  sensors,  and  subsystem 
actuators. 

For  this  reason,  internal  vehicle  three-axis  control  methods  with  lower  mass  are 
needed.  Flywheel-based  systems  providing  both  energy  storage  and  attitude  control 


address  the  need  for  combined  energy  storage  and  attitude  control.  In  particular,  the 
Air  Force  Research  Laboratory  (AFRL)’s  Flywheel  Attitude  Control,  Energy  Trans¬ 
mission  and  Storage  (FACETS)  program  will  combine  all  or  part  of  the  energy  storage, 
attitude  control,  and  power  management  and  distribution  (PMAD)  subsystems  into 
a  single  system,  significantly  decreasing  bus  mass  (and  volume)  by  eliminating  the 
need  for  conventional  chemical  batteries  [2], 

1.1.2  IPACS 

An  Integrated  Power  and  Attitude  Control  System  (IPACS)  such  as  FACETS  employs 
flywheels  as  “mechanical  batteries”  to  perform  the  attitude  control  and  energy  stor¬ 
age  functions.  The  IPACS  concept  eliminates  vehicle  mass  while  improving  system 
performance  and  lifetime.  The  results  of  a  recent  Air  Force  trade  study  are  reflected  in 
figure  1.1.  This  figure  shows  the  45-50  percent  cost  savings  such  a  system  attains  as  it 
removes  the  typical  satellite’s  most  expensive  subsystem,  the  batteries.  Up  until  now, 
the  well-documented  IPACS  concept  was  never  implemented  due  to  high  flywheel  spin 
rates  (on  the  order  of  40K  to  80K  RPM  versus  less  than  5K  RPM  for  conventional 
Control  Moment  Gyroscopes  (CMGs)  or  momentum  wheel  actuators).  At  such  high 
speeds,  the  actuators  quickly  wear  out  mechanical  bearings.  Additional  challenges 
include  flywheel  material  mass/durability  and  stiffness  inadequacies.  Recently,  the 
advent  of  composite  materials  and  magnetic  bearing  technology  has  enabled  realistic 
IPACS  development  [2,  3]. 
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Figure  1.1:  USAF  Trade  Study  Results 


Energy 

Storage 

19% 


Balance  of 
S/C 
45% 


Attitude 

Control 

36% 


Balance  of 
S/C 
71% 


1.1.3  Simultaneous  Control  Difficulty 

The  control  problem  of  simultaneous  energy  storage  and  attitude  control  is  far  from 
trivial,  even  in  its  simplest  conceivable  form.  While  decoupling  the  attitude  control 
and  energy  storage  functions  may  be  a  workable  solution,  research  in  related  areas 
suggests  it  may  not  be  the  best  approach  [2]. 


1.1.4  Momentum  Wheel  Concept 


It  has  been  shown  in  [4,  5]  that  simultaneous  momentum  management  and  power 
tracking  can  be  accomplished  with  four  or  more  wheels  in  reaction  wheel  (RW)  mode. 
This  is  done  by  adjusting  the  wheel  acceleration  in  the  null  subspace  of  the  required 
attitude  control  torque  dynamics  matrix  in  such  a  way  as  to  generate  the  required 
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vehicle  power  while  not  imparting  adverse  external  torque  on  the  spacecraft.  Further¬ 
more,  this  method  was  shown  to  be  practical  for  tracking  the  required  time  history 
profiles  for  several  types  of  satellites  [3,  4]. 

1.1.5  AFRL’s  FACETS 

AFRL’s  FACETS  program  intends  to  employ  VSCMGs  on  the  Advanced  Space  Struc¬ 
tures  Technology  Research  Experiments  (ASTREX)  test  article  (depicted  in  figure 
1.2)  using  a  concept  similar  to  that  for  momentum  wheels  described  in  section  1.1.4. 
This  test  article  rests  on  an  air  bearing  and  provides  one  of  the  best  ground-based 
test  environments  for  simulating  on-orbit  satellite  motion.  Tests  using  this  structure 
produce  results  that  are  very  close  to  that  experienced  by  the  real  space  vehicle. 
The  ASTREX  platform  helps  reduce  satellite  program  costs  since  problems  can  be 
eliminated  long  before  payloads  are  launched  into  orbit. 

Figure  1.3  shows  a  CMG  actuator  configuration  previously  used  on  the  AS¬ 
TREX  structure  in  the  early  1990s.  A  similar  configuration  is  planned  for  the  new 
FACETS  flywheel  actuator  installation.  There  are  different  kinds  of  actuators  that 
could  be  used  for  an  on-orbit  IPACS.  Figure  D.2  shows  a  few  of  these  along  with  each 
alternative’s  associated  specifications. 
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1.2  Problem  Statement 


Simultaneous  attitude  and  power  tracking  for  a  rigid  spacecraft  using  Variable  Speed 
Control  Moment  Gyroscopes  (VSCMGs)  is  possible.  This  work  examines  the  the¬ 
oretical  background  governing  the  use  of  VSCMG  actuators,  presents  a  simulation 
software  tool  designed  to  investigate  systems  using  these  actuators,  and  analyzes  the 
results  of  several  simulation  runs  using  this  software  tool  for  a  simultaneous  attitude 
control  and  energy  storage  system. 

In  this  work,  the  exact  nonlinear  equations,  never  apparently  derived  in  the 
literature,  governing  such  an  attitude  control/energy  storage  system  are  derived  and 
examined.  The  derivation  is  based  on  n  VSCMGs,  n  being  an  arbitrary  number.  The 
resulting  equations  are  independent  of  any  particular  actuator  configuration.  The 
explicit  derivation  enables  direct  application  of  the  theory  to  an  actual  satellite  sys¬ 
tem.  The  generality  of  the  theory  permits  application  to  a  wide  variety  of  spacecraft 
missions.  This  creates  flexibility  for  future  space  systems  contemplating  an  IPACS 
using  VSCMGs  for  attitude  control  and  energy  storage.  In  addition,  the  design  of 
a  SIMULINK-based  software  simulation  tool  for  analyzing  these  types  of  spacecraft 
actuators,  called  the  VSCMG  Workbench,  is  presented.  This  software  will  be  used 
by  AFRL  to  analyze  the  application  of  simultaneous  attitude  control/energy  storage 
to  different  satellite  programs.  Finally,  results  from  several  simulation  runs  using  two 
different  designed  control  algorithms  for  two  different  control  cases  -  attitude  regu¬ 
lation  and  attitude  tracking  -  are  presented.  From  these  controller  scenario  tests  as 
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well  as  overall  system  performance  tests,  the  effects  that  several  parameters  have  on 
the  system  are  analyzed. 

1.3  Research  Outline 

The  structure  of  this  document  supports  the  goals  just  presented  in  section  1.2.  In  or¬ 
der  to  achieve  these  goals,  one  must  investigate  several  aspects  of  the  problem.  These 
aspects  are  decomposed  into  7  chapters  with  further  supporting  details  contained  in 
the  Appendices. 

Chapter  2  provides  a  review  of  the  countless  documents  related  to  this  topic. 
Much  of  this  research  directly  builds  on  the  work  of  others  in  further  developing  the 
theory,  proposing  a  simulation  analysis  tool,  and  testing  the  results.  Therefore,  it 
is  very  important  that  one  is  able  to  put  the  past  work  into  historical  perspective. 
Chapter  2  aims  to  do  just  that. 

Chapter  3  investigates  the  intricate  details  of  the  theory.  This  chapter  is  the 
backbone  upon  which  the  analytical  results  are  based.  It  shows  the  nuts  and  bolts  of 
the  integration  of  previous  results  related  to  spacecraft  attitude  control  and  energy 
storage  using  VSCMGs. 

Chapter  4  presents  the  VSCMG  Workbench,  a  software  simulation  tool  to  ana¬ 
lyze  the  application  of  simultaneous  energy  storage  and  attitude  control  for  satellites. 
The  generality  and  modularity  of  this  tool  renders  it  a  robust  computer  platform  that 
can  be  modified  for  countless  satellite  programs  contemplating  using  the  simultaneous 
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spacecraft  energy  storage  and  attitude  control  methodology  presented. 

Chapter  5  outlines  the  large  battery  of  VSCMG  Workbench  simulation  tests 
performed  to  gain  system  insight.  This  chapter  centers  only  on  the  testing  structure 
and  not  on  the  test  result  analysis. 

Chapter  6  presents  the  results  of  running  the  battery  of  tests  presented  in 
chapter  5.  These  results  illustrate  several  lessons  learned  about  the  example  system 
and  reflect  the  different  parameter  effects  on  the  spacecraft  system  overall  as  well  as 
in  terms  of  its  subsystems.  This  type  of  analysis  allows  users  to  see  what  they  can 
learn  by  running  the  VSCMG  Workbench  and  applying  the  concepts  to  a  particular 
satellite  and/or  satellite  mission. 

Chapter  7  summarizes  the  different  conclusions  one  can  draw  from  this  re¬ 
search.  Thus,  this  chapter  highlights  the  most  important  lessons  learned  during  the 
research  effort. 

Finally,  chapter  8  describes  potential  areas  for  further  study.  This  research 
raises  exciting  questions  about  spacecraft  systems  employing  simultaneous  3-axis  at¬ 
titude  control  and  energy  storage  methods.  Such  further  study  will  help  researchers 
learn  more  about  this  interesting  concept. 
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Figure  1.2:  ASTREX  Test  Article 


Figure  1.3:  CMG  Actuators 
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Attitude  Control 


CMGs  for  2  transverse  axes 


RWA  for  polar  axis 


Honeywell  M325S  CMG 

Nominal  Torque:  440  N-m 
Momentum:  440  N-m-s 
Mass:  65.3  kg 


Honeywell  HR150  RWA 

Nominal  Torque:  1.62  N-m 
Momentum:  203  N-m-s 
Mass:  27.7  kg 


Energy  Storage 

Nickel- Hydro  gen  or  Lithium 
Ion  Battery 

Lithium  lo  n  Battery 

Specific  Energy:  110  W-h/kg 
Depth  of  Discharge:  46% 

Nickel-Hydrogen  Battery 

Specific  Energy:  50  W-h/kg 
Depth  of  Discharge:  56% 


Combined  Attitude  Control  &  Energy  Storage 


FACETS 

N  ominal  T  orque:  580  N-m 
M omentum:  5858  N-m-s 
Specific  Energy:  50  W-h/kg 


Depth  of  Discharge:  82% 


Figure  1.4:  Actuator  Alternatives 
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CHAPTER  2 

LITERATURE  REVIEW 


Combined  spacecraft  attitude  control  and  energy  storage  is  a  well-studied  problem. 
In  fact,  the  concept  of  using  flywheel  based  actuators  for  energy  storage  and  attitude 
control  has  been  around  since  the  1960s  [3].  This  chapter  will  review  the  significant 
literature  related  to  the  problem  of  simultaneous  attitude  control/energy  storage  us¬ 
ing  VSCMGs.  The  relevant  literature  can  be  broken  down  into  two  parts:  i)  flywheels 
used  for  energy  storage  (terrestrially  and  in  space)  and  ii)  control  moment  gyroscopes 
for  attitude  control.  Reference  [3]  gives  an  excellent  overview  of  the  pertinent  litera¬ 
ture  on  the  former  part  and  [6]  remarkably  summarizes  the  latter. 


2.1  Flywheel  Energy  Storage 

The  roots  of  the  concept  of  using  flywheels  to  store  energy  came  to  light  in  the 
1960s  during  investigation  of  high  performance  kinetic  energy  storage  systems  [3]. 
As  Hall  mentions  in  [3], “the  use  of  flywheels  instead  of  batteries  to  store  energy  on 
spacecraft  was  suggested  as  early  as  1961  (by  Roes),  when  a  17  W  hr/kg  composite 


flywheel  spinning  at  10  to  20  thousand  RPM  on  magnetic  bearings  was  proposed. 
The  configuration  included  two  counter-rotating  flywheels,  and  the  author  did  not 
mention  the  possibility  of  using  the  momentum  for  attitude  control.”  Meanwhile,  the 
use  of  reaction  wheels  for  large  spacecraft  precision  pointing  and  attitude  control  has 
been  around  for  well  over  30  years  [7,  8].  As  Hall  states,  once  the  concept  of  using 
flywheels  for  energy  storage  arose,  it  only  made  sense  that  such  flywheels  could  be 
combined  as  momentum  wheels  for  attitude  control  and  energy  storage  [3]. 

Energy  storage  flywheels  rapidly  increased  in  popularity  during  the  late  1970s 
and  early  1980s  as  a  result  of  the  “energy  crisis”  [9].  Several  industries  are  now 
employing  the  use  of  such  flywheels  in  some  fashion  (such  as  for  industrial  plants, 
energy  efficient  automobiles,  and  for  use  in  uninterruptible  Power  Supplies  (UPS))  [3]. 
A  great  conceptual  review  of  flywheel  energy  storage  as  it  applies  to  these  different 
industries  is  given  in  references  [10]  and  [11]. 

Research  related  to  flywheel  energy  storage  and  its  space  applications  con¬ 
tained  in  the  literature  can  be  further  subdivided  into  three  eras:  1970-1980,  1980- 
1990,  and  1990-present.  The  ensuing  sections  will  recount  these  eras. 

2.1.1  Flywheel  Energy  Storage  in  the  1970s 

As  mentioned,  the  concept  of  attitude  control  and  energy  storage  has  been  around 
for  a  while.  Soon  after  the  energy  storage  idea  initially  appeared  in  reference  [12], 
the  combined  concept  arose.  There  were  three  primary  periods  during  which  the 
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combined  concept  was  intensively  addressed.  References  [13,  14,  15,  16,  17]  cover 
the  first  period  of  investigation  (from  1970-1980),  during  which  time  Anderson  and 
Keckler  coined  the  term  IPACS  (Integrated  Power  and  Attitude  Control  System)  to 
describe  such  a  system.  At  that  time,  NASA  did  some  extensive  concept  feasibility 
studies  and  even  held  a  few  working  groups  on  the  topic  in  order  to  investigate 
potential  uses  for  shuttle-age  programs. 

2.1.2  Flywheel  Energy  Storage  in  the  1980s 

References  [18,  19,  20,  21,  22,  23,  24,  25,  26,  27]  encompass  the  second  period  of 
extensive  combined  attitude  control  and  energy  storage  concept  exploration  in  which 
a  Combined  Attitude  Control  and  Energy  Storage  (CARES)  system  was  investigated 
by  the  Charles  Stark  Draper  Laboratory  as  well  as  investigation  into  an  IPACS  for 
the  Space  Station  was  continued  by  NASA. 

The  CARES  system  explored  the  technological  issues  related  to  mounting 
momentum  exchange  devices  on  magnetic  bearings.  Downer,  Eisenhaure,  et  al.  pub¬ 
lished  several  papers  in  the  mid  1980s  that  presented  the  Draper  Laboratory’s  flywheel 
energy  storage  and  spacecraft  attitude  control  system  concept,  magnetic  suspension 
design,  component  specifications,  and  hardware  issues  [19,  20,  22].  This  system  ap¬ 
plied  magnetic  bearing  technology  developed  for  the  Department  of  Energy  to  a  space¬ 
craft  attitude  control  problem  [20].  As  stated  in  [20],  the  resulting  system  provided 
“attitude  control  about  the  roll,  pitch,  and  yaw  axes.”  Reference  [20]  stated  that 
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two  axes  were  “controlled  by  utilizing  the  wheels  as  control  moment  gyros,  while  the 
third  axis  (was)  controlled  by  utilizing  the  wheels  as  reaction  wheels.”  Interestingly, 
the  magnetic  bearing  of  that  system  was  “used  as  a  ‘gimbal’  for  a  control  moment 
gyro”  [20], 

Excellent  summaries  of  the  mid-1980s  technological  issues  facing  IP  ACS  imple¬ 
mentation  were  given  by  Simon  and  Van  Tassel  in  [26]  and  by  Oglevie  and  Eisenhaure 
in  [23].  These  sources  describe  the  material,  durability,  and  stiffness  issues  mentioned 
in  chapter  1  along  with  flywheel  magnetic  suspension  hurdles  to  implementation  and 
the  problem  of  motor,  generator,  and  electronic  technologies  in  the  area  of  energy 
conversion.  The  latter  of  these  references  pointed  out  that  “the  energy  transfer  re¬ 
quired  (for  rapidly  slewing  large  satellites)  is  so  large  that  it  imposes  very  large  peak 
torquing  power  requirements  (many  horsepower).  The  IPACS  capability  to  store  and 
deliver  high  peak  power,  and  regeneratively  brake  the  rapid  slewing  maneuver,  is  ide¬ 
ally  suited  to  these  applications”  [23].  Olmstead  looked  at  applying  flywheel  energy 
storage  technology  for  attitude  control  using  a  counter-rotating  flywheel  concept  [24], 
Flatley  examined  the  use  of  four  flywheels  as  reaction  wheels  in  a  tetrahedron  config¬ 
uration  and  the  associated  voltage  regulation  problem  presented  when  such  a  system 
is  used  for  combined  attitude  control  and  energy  storage  [21].  Studer  and  Rodriguez 
provided  some  of  the  overarching  design  issues  in  implementing  an  Attitude  Control 
and  Energy  Storage  (ACES)  system  [27]. 

As  is  evident,  many  technological  combined  attitude  control  and  energy  stor- 
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age  design  issues  arose  and  were  investigated  in  the  1980s.  Building  on  concepts 
introduced  in  the  early  1970s,  the  research  uncovered  some  of  the  important  design 
barriers  to  IPACS  implementation.  This  paved  the  way  for  re-examination  of  the 
problem  in  the  1990s. 

2.1.3  Flywheel  Energy  Storage  in  the  1990s  and  beyond 

The  third  period  of  study,  addressed  by  references  ([28,  3,  4]),  comprises  recent  in¬ 
vestigation  due  to  the  advance  of  enabling  technologies  mentioned  in  Chapter  1  such 
as  flywheel  durability  and  magnetic  bearing  technology.  The  investigation  found  in 
[4]  includes  a  more  detailed  analysis  of  the  exact,  nonlinear  equations  of  motion  for 
flywheels  used  as  momentum  wheels.  Also,  reference  [29]  is  a  good  source  for  the  Air 
Force  Research  Laboratory’s  current  application  of  flywheel  energy  storage. 

Most  of  the  present  spacecraft  energy  storage  research  has  been  sponsored 
by  two  government  sources,  NASA  (for  future  implementation  on  the  International 
Space  Station)  and  AFRL  in  its  perpetual  effort  to  make  next  generation  spacecraft 
cheaper  and  lighter. 

2.2  Control  Moment  Gyroscopes 

Similar  to  flywheel  energy  storage,  control  moment  gyroscopes  have  been  well-studied 
over  the  last  30  years.  This  research  can  be  broken  into  three  areas:  early  concepts, 
on-orbit  applications,  and  theoretical  analysis.  Each  of  these  areas  will  be  addressed. 
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In  this  research,  we  primarily  concentrate  on  the  use  of  the  Single  Gimbal  CMG 
(SGCMG)  as  this  device  is  much  simpler  than  its  counterpart,  the  Double  Gimbal 
CMG  (DGCMG). 

2.2.1  Early  Control  Moment  Gyroscope  Concepts 

The  concept  of  using  control  moment  gyroscopes  for  attitude  control  arose  in  the  mid- 
1960s  during  the  modern  era  of  space  exploration.  One  of  the  first  sources  that  details 
the  concepts,  benefits,  and  potential  configurations  of  CMGs  for  attitude  control  is 
given  by  Jacot  and  Liska  in  [30].  That  document  reflects  a  commonality  found  in 
several  works  in  the  CMG  literature  -  the  use  of  linearized  equations  of  motion. 
This  work  shows  that  the  beginning  of  the  use  of  CMGs  started  through  laboratory 
investigation  of  the  system  and  stemmed  from  the  fact  that  early  spacecraft  gyroscopic 
sensors  imparted  slight  torques  on  a  spacecraft  during  its  operation  [30] .  Application 
of  these  devices  on  spacecraft  has  mainly  been  for  large  angle,  three-axis  satellite 
maneuvers  due  to  the  CMG’s  mass  and  power  requirements. 

2.2.2  On-Orbit  Uses 

There  have  been  several  large  space  vehicles  that  have  employed  CMG  actuators 
since  the  1960s.  One  spacecraft  that  used  three  orthogonally  mounted  DGCMGs  was 
NASA’s  Skylab  space  station  [8].  Also,  the  Russian  Space  Station  MIR  has  employed 
six  parallel  mounted  SGCMGs  successfully.  Third,  the  International  Space  Station 
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uses  four  parallel  mounted  DGCMGs,  two  of  which  are  mounted  anti-parallel  to  the 
other  two.  But  as  Wie  states  in  [8],  “control  moment  gyros  have  never  been  used  in 
commercial  satellites” ,  perhaps  due  to  the  complex  nature  of  these  systems  and  the 
large  satellites  they  prove  cost  effective  for. 

2.2.3  CMG  Equations  of  Motion  Analysis 

Although  the  use  of  CMGs  for  attitude  control  concept  was  born  in  the  1960s,  it 
was  not  until  the  work  by  Marguilies  and  Aubrun  in  1979  [31]  that  the  first  real 
mathematical  development  of  CMG  theory  arose.  This  work  centered  on  the  related 
geometry  of  using  CMGs  and  introduced  some  of  the  properties  of  CMG  singularities. 

Oh  and  Vadali  [32]  seemed  to  be  the  first  to  seriously  investigate  the  nonlinear 
CMG  equations  of  motion.  References  [33,  34,  35,  36]  further  present  the  singularity 
problems  that  arise  from  the  use  of  Control  Moment  Gyroscopes  in  many  applications. 

2.2.4  Variable  Speed  Control  Moment  Gyroscopes 

The  first  realistic  study  of  combining  the  two  functions  of  reaction  wheel/momentum 
wheels  and  control  moment  gyroscopes  was  presented  by  Ford  and  Hall  in  [6,  37]. 
However,  as  Schaub  points  out  in  [38],  Ford  and  Hall  stop  short  of  investigating  the 
simultaneous  use  of  reaction  wheels  and  control  moment  gyroscopes.  For  this  reason, 
Schaub,  Junkins,  and  Vadali  proposed  the  use  of  Variable  Speed  Control  Moment 
Gyroscopes  in  [38,  39,  40]  to  help  avoid  unfavorable  classical  CMG  singularites  as 


17 


well  as  develop  the  governing  nonlinear  equations  of  motion  associated  with  this  type 
of  system.  The  research  by  Schaub  et.  al.  even  goes  to  the  point  of  analyzing  the 
use  of  null  motion  by  VSCMGs  in  order  to  reorient  a  spacecraft  to  take  advantage  of 
more  robust  preferred  gimbal  angle  configurations. 

2.2.5  The  Present  Research:  VSCMGs  for  Energy  Storage 
and  Attitude  Control 

The  only  apparent  suggestion  of  using  Variable  Speed  Control  Moment  Gyroscopes 
for  combined  attitude  control  and  energy  storage  was  made  in  1985  by  O’Dea  et.  al. 
in  [22]  The  authors  in  that  reference  mentioned  gimballing  the  motor  and  used  the 
term  “variable-speed”  CMGs  while  proposing  that  this  concept  was  “promising”  but 
stopped  short  of  developing  the  governing  the  nonlinear  equations  of  motion.  This 
left  the  door  open  for  research  into  combined  attitude  control  and  energy  storage 
using  Variable-Speed  Control  Moment  Gyroscopes. 

In  summary,  the  present  research  combines  the  previous  results  of  other  in¬ 
vestigations  and  applies  it  via  simulation  to  a  realistic  satellite  system.  Some  of  the 
initial  theoretical  results  of  this  research  are  published  in  [2]  as  well  as  in  an  upcoming 
American  Control  Conference  paper  by  Richie,  Tsiotras,  and  Fausz  [41].  Although 
documentation  in  the  literature  is  very  abundant  for  each  of  these  topics  separately 
(flywheel  energy  storage,  control  moment  gyroscopes,  and  variable-speed  control  mo¬ 
ment  gyroscopes),  the  present  research  combines  them  into  a  more  general  result. 
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The  theoretical  results  of  this  study  are  further  presented  in  the  next  chapter. 
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CHAPTER  3 


THEORETICAL  BACKGROUND 


Next,  the  theory  behind  a  VSCMG-actuated  spacecraft  attitude  control/energy  stor¬ 
age  system  will  be  discussed. 

3.1  Preliminaries 

It  is  necessary  to  present  some  key  definitions  and  mathematical  relations  that  will 
be  useful  in  the  subsequent  sections  so  that  one  can  best  understand  the  theory. 

3.1.1  Notation 

Matrices  are  presented  in  bold  as  in  A,  vectors  are  presented  with  an  arrow  as  in  x, 
and  vector  components  with  respect  to  a  particular  basis  (essentially  used  as  n  x  1 
column  matrices)  are  denoted  in  italics  as  in  x.  The  derivative  of  a  vector  x  with 
respect  to  the  inertial  reference  frame  is  denoted  by 


For  convenience,  in  representing  vector  cross  products,  define  the  skew  symmetric 
matrix  x  £  IR3x3,  for  x  £  IR3,  which  allows  one  represent  x  x  y  =  xy  in  terms  of  a 
column  matrix  of  its  components  in  the  same  basis  as  that  in  which  the  components 
of  are  represented  in  x  and  y.  Thus  the  components  of  x  xy  are  written  as  bx  where 

0  -x3  x2 

x3  0  —x\  (3-2) 

— x2  X\  0 

The  0  operator  represents  the  Kronecker  product  of  two  matrices.  For  any  two 
matrices  A  £  ]Rnxm  and  B  £  JRpxg  the  matrix  A  0  B  is  the  IRm?,><n9  matrix  given  by 

a,\\B  a\2B  •  •  •  a\ nB 

a2\B  a22B  •  •  •  a2nB 

A®  B  = 

&mi  B  am2B  O'mnB 

3.1.2  Assumptions  and  Definitions 

Now,  one  can  make  several  key  assumptions  throughout  the  derivation  of  the  system 
model.  These  assumptions  will  be  identified  individually  in  the  text  as  they  become 
appropriate.  For  the  development  of  the  equations  of  motion,  one  should  consider 
a  system  consisting  of  a  rigid  spacecraft  with  body-fixed  reference  frame,  B,  which 
includes  an  array  of  n  rigid  VSCMGs  with  reference  frames  fixed  to  each  of  the 
VSCMG  gimbals,  Qi,  Q2, . . . ,  Qn.  Figure  3.1  (taken  from  [40])  illustrates  the  Q  frame 
for  one  VSCMG.  As  shown  in  Fig.  3.1,  the  frames  attached  to  the  n  VSCMGs,  Qj, 
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Figure  3.1:  Variable  Speed  Control  Moment  Gyro 


j  =  1, n,  are  characterized  by  an  orthogonal  set  of  unit  vectors,  gSj,  gtj,  and  ggj, 
j  =  1  where  the  subscripts  s,  t,  and  g  denote  the  spin,  transverse  and  gimbal 

axes,  respectively,  satisfying  the  relation 

9gj  x  9s j  —  9tj  (3-3) 

The  matrix  LsGj  €  Hl3x3  is  the  rotation  matrix  transforming  vectors  expressed  as 
components  in  reference  frame  Gj  to  frame  B  and  conversely,  LGjB  =  LJqG]  is  the 
rotation  matrix  from  B  to  Gj.  Conveniently,  the  matrix  LBGj  £  IR3x3  can  be  defined 


as 

Lbgj  =  [fjsj  9tj  9gj ]  (3-4) 

where  the  column  matrices  gsj ,  gtj ,  and  ggj  contain  the  components  of  the  unit  vectors 
9sji  9tj ,  and  gg]  expressed  in  B. 

Next,  let  f]  6  IRn  be  the  column  matrix  that  contains  the  wheel  rotational 
speeds  of  the  n  VSCMGs,  ty,  j  =  1, . . . ,  n,  and  let  7  G  IRn  be  the  column  matrix 
that  contains  the  gimbal  angles  of  the  n  VSCMGs,  7 j  =  1, . . .  ,n.  Then,  0  G  IRn, 
7  G  IRn,  and  7  €  IR”  are  defined  similarly. 

To  simplify  the  presentation  of  the  results,  one  first  defines  Gs  G  IR3*”  as  a 
matrix  whose  columns  are  measure  numbers  of  gSj  in  the  B  basis,  so  that 

Gs  =  [p5i  •  •  •  gsn ]  (3-5) 

and  then  defines  Gt  and  Gg  similarly  for  the  transverse  and  gimbal  axis  unit  vectors, 
respectively.  One  must  also  define  the  matrix  GSd  G  IR3nxn  such  that 

Gsd  =  diag  [gsi,  gS2,  ■  ■  ■ ,  9sn]  (3.6) 

and  similarly  for  Gtd  and  Ggd. 

Next,  it  will  be  convenient  that  one  define  several  matrices  involving  the  inertia 
properties  of  the  VSCMGs.  The  inertia  values  of  each  VSCMG  is  decomposed  into 
the  contributions  of  the  wheel  and  the  gimbal  structure  using  the  scalar  variables  Iwaj , 
Iwtj,  Iwgj,  IgS]  ,  lGtji  and  Icjgj ,  j  =  1, . . . ,  n,  where  the  subscripts  W  and  G  denote  the 
wheel  and  gimbal  structure  contributions  along  the  s,  t ,  and  g  axes,  respectively.  It 
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should  be  noted  that  the  assumption  is  made  (as  is  done  in  [39])  that  all  VSCMGs  are 
perfectly  balanced  and  aligned  so  that  the  unit  vectors  gS] ,  gtj ,  and  gaj,  j  =  1, ...  ,n, 
represent  principal  directions  for  the  VSCMG  reference  frames. 

Now,  the  reader  can  define  the  following  inertia  matrices: 

IGj  =  diag  [lGsj,IGtj,IGgj] 

I Gs  \Jgs\^Gs2  ’  '  ’  ^Gsn J 

IGsi  =  diag  [IGs] 

iGsm  =  lGsd  ®  I3  (3.7) 

where  IGj  G  IR3,  IGs  G  IRlxn,  IGsd  G  IRnx”,  and  IGsm  G  IR3nx3n.  Similarly,  one  can 

define  matrices  for  the  other  gimbal  structure  inertias:  IGt  G  IRlxn,  IGg  G  IRlxn, 
IGtd  G  IRnxn,  IGgd  G  IRnxn,  IGtm  G  IRnxn,  and  IGgm  G  IRnxn;  as  well  as  for  the 
wheel  inertias:  I \y-  G  IR3,  Iws  £  IR*xn,  Iwt  ^  IR^xn,  Iwg  £  IR*xn,  Iwsd  £  IRnxn, 
IWtd  G  Etnxn,  IWgd  G  IRnxn,  IWsm  e  lR3nx3n,  IWtm  e  IR3nx3n,  and  IWgm  G  IR3nx3n; 
and  that  I3  is  the  3x3  identity  matrix. 

At  times  it  is  convenient  to  combine  the  inertia  contributions  of  the  wheel  and 
gimbal  structure,  so  define  Ja  G  nt3x3,  Js  G  IRlxn,  Jsm  <E  IR3nx3n,  and  Jsb  G  IR 3x3,1 
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such  that 


J  j  —  Igj  +  Iwj 
J  s  —  Igs+Iw3 
Jsm  =  lGam  +  Iwsm 
J sb  ~  J s  &1  1 3 

and  similarly  for  Jt  E  IRlxn,  Jg  E  IRlxn,  Jtm  G  IR3nx3n,  Jgm  E  ]R3'ix3'1,  Jtb  e  1R3x3”, 
and  Jgb  E  !R3x3n.  Finally,  one  defines  matrices  Od  E  IRnXn  and  E  IR3nxn  such 
that 


Wrf  =  diag  [u,u, . . .  ,lo]  =  In®  lo  (3.8) 

fld  =  diag  •  •  •  ,^n]  (3.9) 

where  the  spacecraft  body  angular  velocity  vector  to  E  1R3  is  repeated  n  times  in  the 
definition  of  u}d- 

3.2  Dynamics 

In  this  section  the  spacecraft  dynamic  system  model  extending  the  results  of  Oh  and 
Vadali  [32]  to  the  case  of  VSCMGs  is  presented.  More  details  of  this  derivation  can  be 
found  in  Appendix  A.  The  equations  of  motion  are  derived  using  Euler’s  equation  [42] 

j*ys/°  =  ^  ( 'hsys /0)  (3.10) 
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where  hsys/°  is  the  total  angular  momentum  of  the  spacecraft  and  the  VSCMG  cluster 
about  the  combined  system  center  of  mass  (point  O),  given  by 

^sys/0=^P/P’+^P*/0  +  ^W/W* 

+hw*/°  +  Kg/g*  +  hG*/°  (3.11) 

and  is  the  sum  of  the  moments  of  all  external  forces  about  point  0  and  of  all 
external  torques  on  the  spacecraft.  hp/p*  is  the  angular  momentum  of  the  spacecraft 
platform  (without  the  VSCMGs  installed)  about  its  own  center  of  mass  P*.  Similarly, 
hG/G*  is  the  sum  of  the  n-VSCMG  gimbal  structure  angular  momenta  about  each  gim- 
bal  structure’s  center  of  mass,  G*j,  and  hw,/w*  is  the  sum  of  n-VSCMG  wheel  angular 
momenta  about  each  wheel’s  center  of  mass,  Wf,  respectively.  For  the  remainder 
of  this  thesis,  it  is  assumed  that  the  points  G*  and  W*  coincide  for  each  VSCMG. 
Additionally,  hp*/°  is  the  angular  momentum  of  the  center  of  mass  of  the  platform 
(located  at  P*)  with  respect  to  point  O,  hw*/°  is  the  sum  total  of  each  wheel’s  center 
of  mass  angular  momentum  (located  at  point  Gj)  with  respect  to  point  O  ,  and  hGG° 
is  the  sum  total  of  each  gimbal  structure’s  center  of  mass  angular  momentum  (located 
at  point  G*)  with  respect  to  point  O. 

Let  t  be  f sys/,°  expressed  in  B.  Then,  the  system  dynamic  equations  of  motion, 
derived  in  Appendix  A  and  simplified  in  Appendix  B,  may  be  expressed  as  [32] 

t  =  It w  cuPj'W  +  B^j  +  D s' y  +  ECl  +  FG,  (3.12) 

where  u  is  the  angular  velocity  vector  of  the  spacecraft  body  with  respect  to  frame  M 
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expressed  in  frame  B,  r  is  as  defined  above,  and  It  is  the  total  inertia  matrix  given 
by 

N 

It  —  I  sc  +  ^  LBGjJjLTBGj  (3.13) 

j= i 

I sc  is  the  spacecraft  platform  inertia  about  point  O  matrix  plus  the  total  point  mass 
inertia  matrix  of  the  n—  VSCMGs.  The  total  VSCMG  point  mass  inertia  matrix 
can  be  further  subdivided  into  the  total  wheel  point  mass  inertia  matrix  and  the 
total  gimbal  structure  point  mass  inertia  matrix.  The  total  wheel  point  mass  inertia 
matrix  is  the  matrix  obtained  by  summing  the  individual  wheel  point  mass  inertia 
matrices.  Each  wheel  point  mass  inertia  matrix  represents  the  inertia  about  point 
O  of  a  particle  located  at  Wj*  that  contains  the  entire  mass  of  that  wheel.  The 
gimbal  structure  point  mass  inertia  matrix  is  defined  identically  to  that  of  its  wheel 
counterpart.  This  is  expressed  mathematically  as: 

/«  =  /p,0  +  fG''0  +  /w‘/°  (3.14) 

Note  that  the  summation  term  in  (3.13)  is  the  time  varying  portion  of  the  total  inertia 
that  changes  as  the  VSCMGs  move  on  the  spacecraft.  In  addition,  the  coefficient 
matrices  in  (3.12)  are  given  by 


B 

=  Gg  ( IGgd  +  Iwgd) 

(3.15) 

Ds(uj,Q,  7) 

=  D\  +  D2  +  -D3 

(3.16) 

E(  7) 

=  Gs  I Wsd 

(3.17) 

F(u,  7) 

=  U)GS  Iwsd 

(3.18) 
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where 


Di(u,  fi,  7)  -  ( Gt  GTsd  Iwsm  ~  Gs  Gjd  Iwtm )  wd 


+Gt  Iwsd  fid 

(3.19) 

D2(oj) 

-  WGgJgd 

(3.20) 

D3{uj,j) 

=  “  GsGjdIGtm) 

~{GtGTsdJtm  -  GsGTtdJsm)]ud 

(3.21) 

3.3  Kinematics 


One  can  represent  the  spacecraft  orientation  using  Euler  parameters  corresponding 
to  the  transformation  from  the  inertial  reference  frame,  A f,  to  the  vehicle  body  frame, 


B,  as: 


Defining  the  operator  Q(/3)  as: 


Ao 
A 
h 
A 3 


Q(fi)  = 


— Ai  — A2  —A3 
Ao  -  A3  A2 
A3  Ao  — Ai 
-  A2  Ai  Ao 


(3.22) 


(3.23) 
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then  the  spacecraft’s  kinematic  differential  equation  is: 

0=1  «G8)w  (3.24) 

Equations  (3.12)  and  (3.24)  together  allow  for  the  simulation  of  the  spacecraft’s  mo¬ 
tion. 

3.4  Attitude  Stability 

In  this  section  a  control  law  for  attitude  tracking  using  VSCMGs  is  developed.  Note 
that  supplemental  details  of  this  derivation  are  contained  in  Appendix  ??.  One 
begins  by  deriving  a  Lyapunov  based  steering  control  law  as  developed  in  [32]  and 
then  formulates  the  power  tracking  control  in  a  somewhat  analogous  manner  to  that 
done  for  RWs  in  [3,  5]. 

As  in  [32],  a  positive  definite  Lyapunov  function  is  defined  in  terms  of  the 
attitude  errors  /3  —  /5r  and  co  —  ur  as  follows 

V  =  k(/3  —  Pr)T(P  —  Pr)  +  —  ujr)T  It{w  —  0Jr)  (3.25) 

where  fir  and  ojt  are  the  desired  vehicle  reference  attitude  and  reference  angular 
velocity,  respectively.  The  derivative  of  the  Lyapunov  function  can  be  expressed  as 

V  =  —  (w  —  ojr)T  [  kQT(p)Pr  —  It  (<h  —  Cor) 

-  \iT  (u  -  cor)  1  (3.26) 
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It  is  evident  that  V  can  be  made  non-positive  if  one  sets 

kQT(/3)/3r  —  It  (w  —  u>r)  —  \ It  (w  —  ur)  =  K  [ui  —  cor)  (3.27) 

where  JFsT  is  a  positive  definite  gain  matrix.  Since  V  is  non-negative  definite,  the 
resulting  system  is  Lyapunov  stable.  It  can  be  shown  that  V  is  zero  if  and  only  if 
oj  —  ujt  and  j3  =  j3r.  By  LaSalle’s  theorem,  one  can  show  that  the  system  trajectories 
are  stabilized  about  the  desired  reference  attitude.  In  addition,  due  to  the  radial 
unboundedness  of  V,  the  resulting  system  is  actually  globally  asymptotically  stable. 

Next,  notice  that  the  term  \It{u  —  wr)  is  a  linear  function  of  y.  Hence,  one 
can  decompose  this  term  as  follows 

Ry  =  \It{u  —  wr)  (3.28) 

where  the  term  R  is  given  by 

R=\(Jsb~  Jtb)  (GsdGjd,  +  GtdG^d)  (c Crd  ~  OJd)  (3.29) 

If  a  matrix  D  is  defined  such  that 

D  =  DS  +  R  (3.30) 

then  one  can  combine  (3.12),  (3.26),  (3.27),  (3.29),  and  (3.30)  to  yield  the  condition 

By  +  EQ,  +  Dy  +  FQ,  = 

K(lo  —  ur)  —  kQT(P)/3r  —  IroJr  ~  CoItu  +  rd  (3.31) 

where  rd  represents  a  disturbance  torque  on  the  vehicle. 
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3.4.1  Required  Torque  for  Attitude  Tracking 

Next,  if  one  defines  a  column  matrix  containing  the  measure  numbers  of  the  required 
torque  vector,  rr  €  IR3  as 

Tr  —  K(u>  —  ur)  —  kQT(/3)/3r  —  ITujr  —  CjItlo  (3.32) 

then  (3.31)  can  be  written  as 

B'y  -f-  Etl  F£l  -)-  D'y  —  Tr  T  Td  (3.33) 

which  expresses  the  torque  required  for  tracking  in  terms  of  the  physical  parameters 
(or  states)  of  the  system. 

As  shown  in  [39]  if  one  assumes  that  the  gimbal  accelerations  are  small,  then 
the  required  torque  equation  can  be  rearranged  in  terms  of  gimbal  rate  and  wheel  ac¬ 
celeration,  which  represents  the  parameters  typically  controlled  by  commercial  CMGs 
and  RWs,  respectively.  The  resulting  steering  law  (known  as  the  velocity  steering  law) 
is  thus 

EVt  +  Z?7  =  Tr  -  FQ  (3.34) 

3.4.2  Alternative  Control  Law 

An  alternative  control  law  to  that  given  by  equations  (3.33)  and  (3.32)  is  derived 
in  Appendix  .  This  control  law  is  more  intuitive  in  the  case  of  attitude  tracking. 
The  resulting  control  law  is  by  first  defining  the  required  torque  and  then  using  this 
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required  torque  in  equation  (3.33).  The  required  torque  in  this  case  is 

6Pi 

tr  —  Klo  —  KLBRur  —  CoItw  ~  ItLbr^t  +  It^LbroJt  ~  k  §p2  (3.35) 

5(3 8 

where  5(3  is  a  column  matrix  of  errors  in  (3j,j  =  0, 1, 2, 3  and  describes  a  transforma¬ 
tion  between  the  vehicle  body  reference  frame,  B,  and  the  desired  vehicle  reference 
frame  TZ  and  cor  is  the  angular  velocity  of  the  desired  vehicle  orientation  (it  can  be 
thought  of  as  a  virtual  rotating  spacecraft  body  [4])  with  respect  to  inertial  space 
(i.e.  frame  A/"). 

For  the  remainder  of  this  report,  the  control  law  described  in  the  pair  of 
equations  (3.33)  and  (3.32)  will  be  referred  to  as  controller  I;  and,  likewise,  the  control 
law  described  by  equations  (3.33)  and  (3.35)  will  be  termed  controller  II. 

3.4.3  Gimbal  Acceleration  Control 

The  main  advantage  of  a  single-gimbal  CMG  is  its  torque  amplification  property  [30]. 
In  order  to  take  advantage  of  this  property,  one  needs  to  provide  a  velocity  command  7 
to  the  CMG  (and  keep  7  small).  In  fact,  most  standard  CMG  actuators  are  controlled 
via  gimbal  rate  and  not  gimbal  acceleration.  Solving  for  7  directly  from  (3.33)  will 
require  large  gimbal  acceleration  commands  and  hence,  large  gimbal  motor  torques. 
Alternatively,  the  analyst  can  choose  a  velocity  command  7  from  (3.34)  and  then 
implement  this  velocity  steering  law  via  an  outer  control  loop  that  will  keep  the 
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actual  7  close  to  the  desired  gimbal  rates.  This  yields  the  following  equation  for 
the  gimbal  acceleration  command  (assuming  the  desired  gimbal  acceleration,  7^,  is 
negligible) : 

7  =  ^  {id,  ~  7) ,  A  >  0  (3.36) 


3.5  Power  Tracking 


The  kinetic  energy  7)  of  the  jth  actuator  is  given  by: 


Twj  =  loJwjlwj^wj 


(3.37) 


where  u>wj  represents  the  angular  velocity  of  the  wheels  with  respect  to  the  gimbal 
structure  written  in  Qj  components.  This  can  be  re-written  as: 


U>wj  — 


0 


(3.38) 


Next,  note  that  for  n  actuators,  the  total  energy  is  just  the  sum  of  each  of  the 
individual  actuator  energies 

n 

Tw  =  YJTWj  (3.39) 

j'=i 

As  alluded  to  earlier,  taking  the  first  derivative  of  the  energy  yields  the  power 
generated  by  the  wheels 


P\v  =  I  wsd&  (3.40) 
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3.6  Simultaneous  Attitude  and  Power  Tracking 


The  available  control  inputs  are  the  rate  of  change  of  the  wheel  speeds  and  the  angular 
velocity  of  the  gimbals  (for  the  case  of  gimbal  velocity  steering  law).  Hence, 

f i  =  «mw  (3.41) 

'j  ~  ^cmg  (3.42) 


where  umw  £  IRn  and  ucmg  £  IRn  are  the  control  inputs  in  the  “momentum  wheel” 
and  “CMG”  modes,  respectively.  Let  the  combined  control  input  u  =  [it^w  uJmg]T. 
Then  after  some  algebraic  manipulation  one  can  write 


where 


“ 

- 

Cn 

C12 

U  = 

F 

x  pv 

Cn 

c22 

1 _ 

(3.43) 


C  u  —  E,  C 12  —  D 

(3.44) 

=  nTiwsd>  C22  —  0  n 

(3.45) 

Furthermore, 

Fpv  =  tt  —  FFl  and  Pfv  =  Pw  (3.46) 


where  Fpv  is  as  defined  in  (3.34)  and  Pfv  is  defined  in  (3.40).  This  system  of  equations 
is  similar  to  the  velocity  steering  laws  defined  in  the  literature  [32,  40],  which  involve 
solving  for  7. 
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Defining  C\  =  [Cu  C 12]  and  C2  =  [C21  C22],  the  equation  for  from 


(3.43)  can  be  written  as 


C\U  —  Fp 


The  general  solution  to  (3.47)  is  given  by 


Vi  -  C-yPjV  T  Vn 


(3.47) 


(3.48) 


where  the  symbol  f  denotes  the  suitable  generalized  inverse,  and  C\Un  =  0  (i.e.,  un 
is  in  the  null  space  of  C 1,  N(Ci)).  Now  substitute  (3.48)  into  the  equation  for  Pfv 
from  (3.43)  so  that 

C2U  =  C2(C\Fpv  +  un )  =  Pfv  (3.49) 

and 

C2un  =  Pm  (3.50) 

where  Pm  —  Pfv  —  C2C\Fpv.  Since  un  E  J\f(Ci),  one  can  find  a  vector  u  such  that 

Un  =  P  nv  (3.51) 

where  P n  =  In  —  C\C\  is  the  orthogonal  projection  onto  J\f(C  1).  Then,  from  (3.50) 
and  (3.51),  and  making  use  of  the  fact  that  Pn  is  a  projection  matrix,  one  can  see 
un  is  chosen  such  that 

un  =  PNCl{C2PNCl)-lPm  (3.52) 

This  completes  the  solution  for  u  of  (3.48)  for  combined  attitude  and  power  tracking. 

In  summary,  given  the  reference  attitude  to  track,  /9r  and  oor,  the  required 
power  Pw  and  the  state  of  the  system  /3,  u>,  Q  and  7,  one  calculates  the  required 
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attitude  tracking  torque  rr  from  (3.32)  and  the  corresponding  control  inputs  umw 
and  iiCgm  from  (3.48)  and  (3.52).  If  a  gimbal  acceleration  steering  law  is  required 
(to  command  the  gimbal  motors),  then  equation  (3.36)  must  be  used  to  “back-step” 
the  velocity  command  «cmg  =  %  to  an  acceleration  command.  Interestingly,  the 
formulated  power  equation  does  not  change  regardless  of  these  steering  laws  as  7  and 
7  do  not  appear  in  the  power  equation. 


36 


37 


CHAPTER  4 

VSCMG  WORKBENCH 
LAYOUT  AND  DESIGN 


The  VSCMG  workbench  is  a  software  program  designed  for  studying  the  control 
system  design  of  a  spacecraft  equipped  with  VSCMGs  for  combined  attitude  control 
and  energy  storage.  Its  modularized  structure  permits  flexibility  in  implementation. 

The  several  layers  of  the  design  are  described  here  as  well  as  basic  program 
operation.  Note  that  a  couple  more  detailed  schematics  of  the  system  design  can  be 
found  in  Appendix  F. 


4.1  Main  Program  Control  Module 

The  first  block  the  user  needs  is  the  Main  Program  Control  Module  shown  in  fig¬ 
ure  4.1.  In  this  module,  the  user  loads  the  reference  data  (which  currently  involves 
tailoring  a  matlab  *.m  file  where  the  user  defines  2  column  matrices  (one  contains  the 
time  data,  the  other  contains  the  parameter  time  history)  for  each  of  the  11  reference 


values  (these  values  are  the  4-parameter  Euler  parameter  history,  3-parameter  vehicle 
angular  velocity  (given  in  the  vehicle  coordinate  reference  frame,  commonly  referred 
to  as  p,  q,  and  r),  and  the  3  angular  velocity  derivatives  (typically  referred  to  as  p ,  q, 
and  f)). 


whole_sim 


Figure  4.1:  Main  Program  Control  Module 
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4.2  Feedback  System  Overview 


Underneath  the  Program  Control  Module  mask  resides  the  overall  feedback  control 
system.  Unless  the  user  intentionally  reveals  the  support  code  by  selecting  to  look 
under  the  SIMULINK  mask  all  underlying  subsystems  are  masked  from  the  user.  The 
underlying  system  includes  modules  representing  characteristics  inherent  in  a  typical 
control  system  -  the  plant  module,  the  plant  actuators,  the  controller,  the  sensor 
suite,  and  the  reference  values.  This  is  reflected  in  figure  4.2.  Each  one  of  these 
major  subsystems  is  introduced  in  this  chapter. 


|)Lroeas|« - 

meas  ure  ment_scopes 


x_ext_meas 

JLintjneas 

xact 

*_sen$e_p!ot 

Sensors 

Figure  4.2:  System  Schematic 
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4.3  System  Parameters  Module 


The  first  underlying  subsystem  is  the  System  Parameters  Module  reflected  in  figure 
4.3.  There  are  two  primary  sets  of  parameters  that  drive  the  system  configuration. 
These  are  set  by  the  user  and  depend  on  the  number  of  actuators  used.  These 
parameters  include  the  reference  VSCMG  Q  frame  installation  orientation  and  the 
system  inertias  (as  described  in  section  3.2). 


4.4  Reference  Module 

Besides  setting  the  system  parameters,  another  key  set  of  inputs  to  the  control  system 
is  the  reference  attitude  and  the  reference  power  that  the  control  system  is  tasked  to 
achieve.  The  flexibility  in  the  design  allows  different  reference  attitude  and  reference 
power  time  histories  (in  the  form  of  data  tables)  to  be  used  as  input  to  the  system. 
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These  can  represent  vehicles  in  different  configurations,  different  vehicles,  different 


actuator  systems,  and  the  list  goes  on. 


Reference  Attitude 


Figure  4.4:  Attitude  and  Power  Reference 


4.5  Controller 

The  central  “brain”  of  the  control  system  is  the  controller.  The  crux  of  the  developed 
theory  directly  relates  to  the  logic  employed  in  the  controller  (see  4.5  and  4.6).  The 
controller  uses  the  vehicle  configuration  mentioned  in  the  previous  sections  along 
with  the  desired  attitude  and  power  profiles  to  determine  stable  actuator  torques 
that  will  steer  the  vehicle  to  the  desired  attitude  while  simultaneously  meeting  the 
power  tracking  requirements  and  heeding  the  torque  amplification  constraint. 
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4.6  Actuators  Module 

The  commanded  vehicle  torque  is  generated  by  the  actuators  shown  in  figure  4.7. 
The  actuators  are  divided  into  two  major  portions:  the  wheel  motors  and  the  gimbal 
motors.  Both  systems  (the  wheels  and  the  gimbals)  turn  commanded  motor  torque 
voltages  into  actuated  motor  torques.  These  torques  are  used  to  calculate  the  actuated 
torque  about  the  combined  system  center  of  mass  (point  O  as  described  in  Chapter 
3).  As  in  several  of  the  other  modules,  this  module  has  been  designed  to  allow  future 
designers  to  add  model  fidelity.  This  way,  more  realistic  models  of  the  actuator  motors 
can  be  added  to  the  simulation  in  order  to  analyze  different  actuator  behaviors  and 
concepts  of  control. 

4.7  Plant  Module 

The  plant  houses  the  modeled  vehicle  inertias,  the  disturbance  model  inputs,  the 
actuated  torque  computation  function,  and  the  integrated  vehicle  state  equations 
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Figure  4.6:  Controller  Parameters  Overview 
(i.e.  for  vehicle  rate  and  vehicle  acceleration  given  the  torques  acting  on  the  vehicle) . 


4.8  Sensors  Module 

Like  the  actuators  module,  the  sensor  module  allows  much  room  for  the  user  to 
add  more  realistic,  higher-fidelity  models.  Such  models  may  include  typical  sensor 
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Figure  4.7:  Actuator  Motor  Transfer  Functions 


Figure  4.8:  Plant  Overview 


inefficiencies,  sensor  data  non-linear  sensor  phenomena  (such  as  saturation,  hysteresis, 
and  backlash),  and  sensor  fusion/weighting  algorithms.  All  these  topics  permit  much 
further  in-depth  analysis  of  on-orbit  spacecraft  phenomena  never  before  feasible. 


4.9  Summary  of  Program  Operation 

In  order  to  use  these  different  modules,  the  must  have  a  good  idea  of  how  the  simu¬ 
lation  works.  As  far  as  basic  program  operations,  the  user 

1.  Determines  where  the  actuators  are  installed  (i.e.  where  the  initial  gimbal 
reference  axis  for  each  VSCMG  is  located  relative  to  the  body  reference  frame). 
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Figure  4.9:  External  System  State  Overview 

2.  Sets  the  vehicle  platform  plus  actuator  point-mass  inertia  matrix  as  well  as 
the  individual  actuator  inertia  matrices  (with  respect  to  each  Qi  frame).  These 
values  are  set  in  the  system  parameters  module. 

3.  Configures  the  controller  parameters  (as  appropriate)  in  the  controller  subsys¬ 
tem  module. 


Figure  4.10:  Sensor  System 
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4.  Sets  the  desired  actuator  model  in  the  actuators  system  module  according  to 


the  desired  fidelity/response  of  the  actuated  system. 

5.  Runs  the  load  ref  function  which  will  load  reference  data  according  to  the 
desired  scenario. 

6.  Starts  the  simulation  in  the  Master  program  module  by  selecting  the  play 
button. 

7.  Selects  generate  plots  in  order  to  create  MATLAB  output  plots,  if  desired, 
when  the  simulation  is  done.  Note  that  theload  ref  and  generate  plots 
functions  should  be  tailored  according  to  the  purpose  to  which  this  program  is 
applied. 
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Figure  4.11:  Disturbance  Torque  Inputs 
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CHAPTER  5 

SIMULATION  TESTS 


Several  different  simulation  tests  were  performed  in  order  to  study  simultaneous  space¬ 
craft  attitude  control  and  power  tracking.  These  tests  can  be  broken  down  into  several 
key  areas:  controller  comparison  tests,  controller  parameter  tests,  system  parameter 
tests,  simulation  parameter  tests,  and  actuator  model  tests. 

Before  describing  these  tests  in  more  detail,  one  must  first  review  some  pre¬ 
liminary  information  on  control  mode  weighting  and  singularity  avoidance  and  then 
outline  the  scenarios  used  for  attitude  regulation  and  attitude  tracking.  Note  that  in 
this  chapter,  only  the  test  scenarios  and  parameters  are  presented  -  analysis  of  the 
tests’  results  follows  in  the  ensuing  chapter.  Furthermore,  it  is  important  to  under¬ 
stand  that  the  baseline  tests  described  for  each  test  group  essentially  serve  as  the 
“experimental  control”  variables  and  thus  allows  one  to  compare  different  parameter 
changes  with  the  baseline  in  order  to  evaluate  a  parameter’s  impact  on  the  system 


performance. 


5.1  Weighting  and  Singularity  Avoidance 


Define  the  weighted  generalized  inverse  C\  identical  to  that  used  in  [40] 

C[  =  WC\T  (CiWCxT)~l  (5.1) 

where  W  is  a  diagonal  RW/CMG  mode  weighting  matrix  such  that  Wgj0  represents 
weighting  the  predominant  CMG  mode  (to  capitalize  on  its  torque  amplification  prop¬ 
erty)  and  is  constant  throughout  a  maneuver.  WSj  is  the  RW  mode  weight  that  comes 
into  play  near  a  CMG  singularity.  WSj  is  given  by 

Wsj  =  Wsj  oexp"^  (5.2) 

in  which  /j  and  Wsjo  are  chosen  by  the  control  system  designer  to  elicit  the  desired 
performance.  5  is  a  parameter  that  describes  the  proximity  to  a  CMG  singularity 
[37,  32,  39,  40].  The  expression  for  5  used  here  is  different  from  that  of  [40]  and  is 
defined  as  the  minimum  singular  value  of  C\ .  This  is  a  more  accurate  way  to  describe 
the  singularity  of  the  matrix  C\  [43]. 


5.2  Attitude  Regulation  Scenario 

The  attitude  regulation  scenario  is  based  on  the  examples  found  in  [32],  [40],  and  [39]. 
Parameters  for  the  baseline  attitude  regulation  scenario,  which  uses  controller  I  from 
Chapter  3,  are  included  in  table  5.1 

Similarly  to  these  references,  a  standard  four  VSCMG  pyramid  configuration 
is  used  here,  which  is  one  in  which  the  VSCMGs  are  installed  so  that  the  four  gimbal 
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Table  5.1:  Regulation  Scenario  Parameter  Settings 


Symbol 

Value 

Units 

N 

4 

unitless 

e 

54.75 

deg 

■son 

[0  0  0] 

rad/sec 

HI 

[0.5  0.5  0.5  0.5] 

unitless 

LOr 

[0  0  0] 

rad/sec 

Pr 

[1  0  0  0] 

unitless 

f  —  _2L  _2L  21  1 

U  4  4  4  J 

rad 

HU 

[0  0  0  0] 

rad 

H 

[50000  60000  55000  65000] 

RPM 

wsj0 

40 

unitless 

Wgj  o 

1 

unitless 

IWj 

diag  {0.70,  0.20,  0.20} 

Kg  m2 

1 Gj 

diag  {0.10  0.10  ,0.10} 

Kg  m2 

I  SC 

diag  {15053,  6510,  11122} 

Kg  m2 

K 

diag  {700,  700,  700} 

Kg  m2/sec 

k 

35.0 

Kg  m2/sec 

k 

10-4 

unitless 

A 

1 

unitless 
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axes  form  a  pyramid  with  respect  to  the  body.  The  pyramid  configuration  has  been 
implemented  here  in  order  to  facilitate  comparisons  with  the  related  literature  even 
though  the  theory  applies  generically  to  the  n-actuator  case.  Figure  5.1  which  was 
taken  from  reference  [38]  depicts  a  standard  CMG/VSCMG  pyramid  configuration. 
Note  that  6  represents  the  pyramid  angle  of  each  VSCMG  that  is  measured  from  the 


Figure  5.1:  VSCMG  Pyramid  Configuration 

vehicle’s  b\  —  62  plane  to  the  VSCMG’s  gimbal  axis. 

In  the  regulation  scenario,  the  goal  is  to  bring  the  vehicle  from  some  arbitrary 
attitude  to  the  reference  attitude  while  tracking  a  given  power  profile.  The  reference 
attitude  in  this  case  is  assumed  to  be  a  fixed  orientation  with  respect  to  frame  M. 
Without  loss  of  generality,  one  can  assume  this  orientation  (in  the  regulation  case)  is 
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the  origin  as  is  done  in  [32]  and  [40]  since  one  can  always  redefine  the  origin  as  being 
the  desired  inertial  orientation. 


5.3  Attitude  Tracking  Scenario 

The  tracking  scenario  is  based  on  the  example  used  in  [4],  which  is  briefly  reviewed  in 
Appendix  E.  Similar  to  the  attitude  regulation  scenario  presented  in  section  5.1,  the 
goal  in  the  attitude  tracking  case  is  to  bring  the  vehicle  from  an  undesired  attitude 
to  the  reference  attitude  while  tracking  a  given  power  profile.  The  given  reference 
attitude  position,  velocity,  and  acceleration  correspond  to  an  on-orbit  example  similar 
to  that  used  in  [4]  in  which  a  near-polar  orbital  satellite  has  to  meet  specific  sun  and 
ground  tracking  requirements. 

The  parameters  used  in  this  scenario  are  virtually  identical  to  that  used  for  the 
regulation  case  (see  table  5.1).  The  only  difference  is  that  instead  of  regulating  about  a 
fixed  angular  position,  angular  velocity,  and  angular  acceleration,  the  Iridium  25578 
orbital  parameters  are  used  to  generate  time  varying  reference  values  for  angular 
position,  angular  velocity,  and  angular  acceleration.  The  relevant  parameters  are 
summarized  in  table  5.2. 
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Table  5.2:  Tracking  Scenario  Parameters  Settings 


Value 

Units 

N 

4 

unitless 

e 

54.75 

deg 

w(0) 

[0  0  0] 

rad/sec 

m 

[0.5  0.5  0.5  0.5] 

unitless 

U)r 

Iridium  25578 

rad/sec 

fir 

Iridium  25578 

1"  7C  7T  77  77  1 

L4  “4  “4  4  J 

rad 

7(0) 

[0  0  0  0] 

rad 

0(0) 

[50000  60000  55000  65000] 

RPM 

wsj  o 

40 

unitless 

Wgj  o 

1 

unitless 

IWj 

diag  {0.70,  0.20,  0.20} 

Kg  m2 

IGj 

diag  {0.10  0.10  ,0.10} 

Kg  m2 

I  SC 

diag  {15053,  6510,  11122} 

Kg  m2 

K 

diag  {700,  700,  700} 

Kg  m2/sec 

k 

35.0 

Kg  m2/sec 

V 

10"4 

unitless 

A 

1 

unitless 
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5.4  Significant  Tests 


The  first  set  of  tests  compares  the  responses  of  the  two  controllers  presented  in  chapter 
3  for  both  scenarios,  regulation  and  tracking.  From  it,  one  can  see  how  effective  each 
controller  is  for  each  scenario. 

Next,  there  are  a  few  controller  parameters  one  has  the  freedom  to  select  when 
using  one  of  the  presented  control  laws.  This  set  of  tests  outlines  the  effects  that  these 
parameters,  the  matrix  gain,  K,  the  scalar  gain  k,  the  reaction  wheel  mode  weighting, 
Wso,  and  the  singularity  scaling  variable,  ju,  have  on  the  response  of  the  system. 

Third,  in  order  to  best  understand  the  spacecraft  system,  one  needs  a  grasp 
of  how  the  different  system  variables  affect  the  system  response.  The  primary  system 
parameters  tested  here  are  the  vehicle  body  and  actuator  inertias.  These  inertias 
have  differing  affects  on  the  controlled  system  response.  One  might  note  that  another 
big  factor  affecting  the  system  is  the  location  of  the  actuators  (i.e.  the  actuator 
configuration) .  A  study  of  the  actuator  locations  is  beyond  the  scope  of  this  research 
but  the  possibilities  are  endless. 

Fourth,  the  simulation  of  a  large  space  structure  of  this  magnitude  and  com¬ 
plexity  obviously  has  several  parameters  related  to  the  simulation  design  and  perfor¬ 
mance  that  impact  the  realism  of  the  simulation.  In  this  set  of  tests,  one  is  able  to 
get  a  feel  for  how  simulation  parameters  such  as  integration  method,  time  step,  and 
run  time  affect  the  realism  of  the  simulation. 

Finally,  tests  were  performed  to  examine  the  effects  of  modeling  the  actuators. 
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The  actuators  were  modeled  in  both  as  identity  gain  and  as  non-identity  gain  models. 
As  in  the  other  tests,  the  results  are  presented  in  chapter  6. 

5.5  Other  Tests  Considered 

Several  other  tests  were  run  during  this  research  effort.  However,  many  of  them  did 
not  appreciably  affect  the  outcome  of  the  simulations  and  have  thus  been  left  out  of 
this  analysis. 

This  includes  issues  such  as  disturbance  torque  effects,  non-diagonal  actuator 
and  plant  inertia  matrices,  and  gimbal  structure  inertia  changes.  The  latter  of  these 
was  tested  but  there  were  very  minor  affects  on  the  overall  system.  The  disturbance 
torques  were  only  initially  tested  and  a  full  study  of  system  disturbance  rejection 
could  be  done  in  the  future.  Third,  incorporating  non-diagonal  inertias  was  not  fully 
modeled  in  the  plant.  It,  too,  could  be  examined  in  future  research. 

5.6  Test  Chapter  Summary 

This  chapter  introduced  several  simulation  tests  carried  out  during  this  research. 
Although  intended  to  cover  a  large  portion  of  the  satellite  system  envelope,  this  set 
of  tests  is  by  no  means  a  completely  exhaustive  suite  of  tests.  Its  intention  is  to  give 
the  satellite  control  system  designer  more  examples  of  different  system  effects.  The 
results  analysis  of  this  battery  of  tests  is  presented  in  chapter  6. 
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CHAPTER  6 

ANALYSIS  OF  TEST  RESULTS 


Once  the  simulation  tests  run  for  this  research  are  understood,  one  can  examine  the 
test  results.  As  was  true  of  the  test  descriptions  (see  chapter  5) ,  the  test  results  can 
be  broken  down  into  several  key  areas:  controller  comparison  results,  control  param¬ 
eter  effects,  system  parameter  effects,  simulation  parameter  effects,  and  actuator  test 
results. 


6.1  Controller  Comparison  Results 

In  this  section,  one  is  able  to  see  the  basic  simulation  output  of  the  scenarios  described 
in  chapter  5  and  the  scenario  data  tables,  tables  5.1  and  5.2. 

6.1.1  Regulation 

The  first  series  of  plots  illustrates  the  regulation  scenario  output.  Figures  6.1  and 
6.2  show  the  results  of  using  controller  I  and  controller  II,  respectively,  which  were 


introduced  in  chapter  3. 


These  plots  (from  left  to  right  starting  at  the  top)  show  time  histories  of  the 
vehicle  angular  position  (in  terms  of  the  Euler  parameter  representation),  angular 
position  error,  angular  velocity,  angular  velocity  error,  the  actuator  wheel  speeds, 
actuator  gimbal  angles,  actuator  gimbal  rates,  and  the  singularity  parameter  for  con¬ 
troller  I.  A  similar  set  of  plots  shows  the  time  response  of  controller  II. 

Figure  6.1  provides  a  good  demonstration  of  the  typical  system  behavior  based 
on  the  given  initial  conditions.  This  plot  shows  the  system  going  to  the  desired 
attitude  and  regulating  motion  about  it  while  simultaneously  tracking  the  given  power 
profile  (shown  in  figure  6.5).  One  can  see  from  these  baseline  plots  the  similarities  in 
the  controllers. 

However,  figures  6.3  and  6.4  show  a  more  direct  comparison  of  the  controllers 
in  terms  of  the  different  components  of  the  angular  velocity  error,  angular  position 
error,  and  the  control  inputs.  Note  that  this  channel-by-channel  plotting  format  is 
used  extensively  throughout  the  remainder  of  this  chapter  in  order  to  present  the 
effects  that  multiple  parameters  have  on  the  system  output.  The  results  are  further 
broken  down  by  vehicle  response  (i.e.  the  components  of  the  angular  velocity  error, 
co  —  u)r  for  controller  I  and  5co  for  controller  II  -  and  the  angular  position  error  -  f3  —  j3r 
for  controller  I  and  5/3  for  controller  II)  This  presentation  method  permits  helpful 
insight  into  the  effects  that  certain  parameters  have  on  the  overall  system  over  time. 

Figures  6.1  and  6.2  show  that  controller  I  and  controller  II  are  very  similar  in 
terms  of  speed  of  response.  The  primary  difference  appears  to  be  that  the  stabilizing 
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vehicle  trajectory  path  is  different  for  each  controller. 

In  figures  6.3  and  6.4,  the  blue  line  for  each  plot  represents  the  controller  I 
time  response  for  each  channel.  The  first  column  of  sub  plots  in  figure  6.3  contains  the 
angular  velocity  sub-component  error  time  histories  and  the  second  column  contains 
the  angular  position  error  time  response.  In  figure  6.4,  the  actuator  responses  given 
for  each  controller  scenario  are  further  broken  down  into  three  columns  corresponding 
to  the  respective  wheel  speed,  gimbal  angle,  and  gimbal  rate  histories  of  each  VSCMG. 
In  the  cases  of  both  plots,  the  blue  line  is  the  controller  I  response  and  the  green  line 
is  the  controller  II  response. 

Thus,  it  is  clear  from  this  plot  that  although  the  response  performance  is 
comparable  between  the  two  controllers,  in  terms  of  vehicle  response,  controller  II 
provides  a  smoother,  more  even  result.  Likewise,  in  comparing  the  actuator  responses 
displayed  in  6.4,  it  is  evident  that  the  only  difference  in  the  resulting  vehicle  motion  is 
reflected  in  the  gimbal  angles.  It  appears  in  this  plot  that  controller  II  provides  a  less 
oscillatory  result.  Since  highly  oscillatory  motion  usually  has  negative  impacts  on  a 
structure’s  lifetime  and  effectiveness,  for  this  scenario  controller  II  is  more  favorable 
than  controller  I. 

In  addition,  one  can  see  in  figure  6.5  the  power  tracking  response  based  on 
the  power  profile  given  in  section  5.2.  As  the  power  tracking  response  is  identical  for 
both  controllers,  only  one  actual  power  output  versus  required  power  profile  plot  is 
included. 
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Note  that  it  is  difficult  to  distinguish  the  actual  response  from  the  desired 
profile.  This  makes  sense  as  for  these  test  runs,  it  was  assumed  that  the  actuator 
motors  produce  exactly  the  desired  motor  output  torque.  Obviously,  as  one  adds 
fidelity  to  the  actuator  models  and  these  models  become  more  realistic,  there  will 
be  some  minor  differences  in  the  desired  and  output  power  as  no  actuator  system  is 
completely  perfect.  Actuator  modeling  affects  are  mentioned  later  in  this  chapter. 

6.1.2  Tracking 

One  can  see  the  baseline  tracking  scenario  output  histories  for  the  two  controllers 
in  figures  6.6  and  6.7  based  on  the  tracking  scenario  mentioned  in  chapter  5  and 
appendix  E.  The  same  power  profile  as  that  used  for  regulation  is  also  used  here. 

The  plots  in  figures  6.6  and  6.7  help  one  gain  a  feel  for  how  the  controllers 
perform  as  far  as  attitude  tracking  is  concerned.  As  before,  both  controllers  track 
the  desired  attitude  in  about  the  same  amount  of  time.  However,  these  plots  indicate 
that  controller  II  provides  a  slightly  more  damped  response. 

The  next  figures  (6.8  and  6.9)  show  the  direct  comparison  of  these  two  con¬ 
trollers.  The  results  of  this  scenario  are  very  much  like  those  of  the  regulation  scenario 
as  is  amplified  in  these  plots.  In  other  words,  the  tracking  scenario  seems  to  bring 
out  the  smoother  response  of  controller  II  over  controller  I.  Hence,  for  similar  reasons 
as  before,  it  is  clear  that  controller  II  is  better  in  this  situation. 
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6.2  Control  Parameter  Effects 


Next,  the  effects  that  the  controller  parameters,  i.e.  the  matrix  gain  K,  the  scalar 
gain,  k,  the  momentum  wheel/control  moment  gyroscope  weighting,  Ws0,  and  the 
singularity  parameter  scaling  constant,  fj,,  have  on  the  system  can  be  examined. 

Figure  6.10  includes  the  vehicle  response  to  changes  in  the  mentioned  control 
parameters.  In  this  figure,  the  dark  green  line  in  the  subfigures  represents  the  vehicle 
response  with  a  large  change  in  the  scalar  gain,  k.  Figure  6.11  reflects  the  impact 
perturbations  in  the  control  parameters  have  on  the  actuators. 

In  these  two  figures  (6.10  and  6.11),  the  forest  green  line  represents  a  large 
change  in  the  scalar  gain,  k.  The  magenta  line  is  the  system  response  to  a  large  change 
in  the  norm  of  K.  The  red  line  relates  to  a  vehicle  simulation  with  a  large  decrease 
in  the  singularity  parameter  scaling  constant  /i  (which  the  control  designer  selects 
based  on  the  desired  performance),  the  light  blue  line  represents  an  increase  in  the 
norm  of  the  momentum  wheel  spin  weighting  WS0)  and  the  blue  line  is  the  regulation 
controller  I  comparison  baseline.  Here,  the  blue,  light  blue,  and  red  lines  in  the  sub 
figures  are  extremely  close  and  therefore  difficult  to  discriminate.  On  the  other  hand, 
the  forest  green  line  is  clearly  much  more  oscillatory  than  the  other  responses.  This 
yields  the  slightly  exaggerated  but  none  the  less  instructive  effect  of  increasing  k.  It 
was  increased  by  an  order  of  magnitude,  from  a  value  of  35  to  350  in  this  example. 
On  the  other  hand,  increasing  the  norm  of  K  resulted  in  a  more  damped  response 
that  is  centered  about  the  baseline  response  (i.e.  the  blue  line).  This  shows  that  as 
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this  gain  matrix’s  norm  is  increased,  the  result  is  a  more  damped  response.  Thus,  as 
is  common  in  a  feedback  control  system,  one  can  manipulate  the  system  performance 
speed  and  damping  through  selection  of  the  feedback  gains  (k  and  K). 

6.3  System  Parameter  Effects 

Similar  to  studying  the  control  parameters,  one  can  analyze  the  impact  that  system 
parameters  (i.e.  system  inertias)  have  on  the  resulting  vehicle  motion.  These  effects 
are  given  in  figures  6.12  and  6.13. 

In  these  plots,  the  blue  line  is  again  the  baseline  controller  I  during  regulation 
response  and  the  forest  green  line  the  vehicle  response  due  to  an  increase  in  vehicle 
platform  inertia  plus  actuator  point-mass  inertia.  The  red  line  represents  increasing 
the  norm  of  all  the  wheel  inertia  matrices;  the  light  blue  line  represents  increasing  all 
the  wheel  spin  axis  inertia  values;  and  the  magenta  line  reflects  increasing  only  one  of 
the  wheel  spin  axis  inertias.  It  should  be  noted  here  that  three  of  the  five  compared 
runs  were  only  simulated  from  0  to  500  seconds,  so  this  is  the  region  of  primary  interest 
on  these  plots.  Here,  one  can  see  a  large  spike  in  the  magenta  response  where  that 
wheel’s  motion  was  adjusted  due  to  the  larger  inertia.  These  results  show  the  vehicle 
going  to  the  target  state  but  at  the  cost  of  some  excessive  wheel  motion  as  compared 
to  the  baseline  case.  It  should  be  noted,  though,  that  this  result  mainly  occurs  when 
modifying  any  and  all  of  the  wheel  spin-axis  inertia.  As  these  inertia  values  are  much 
higher  than  the  other  axis  inertia  or  the  gimbal  structure  inertia  (this  is  clearly  due 
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to  the  assumed  circular  disk  wheel  shape  for  the  flywheels).  One  can  also  see  in  the 
actuator  plots  that  wheel  spin-axis  inertia  changes  result  in  more  oscillatory  actuator 
(and  thus  vehicle)  motion  as  compared  to  the  baseline.  Interestingly,  as  reflected  in 
the  red  line’s  performance,  if  all  the  wheel  inertias  are  changed  (e.g.  scaled  sizing 
along  all  axes  of  the  wheel  system),  then  the  resulting  motion  is  more  favorable.  Also 
interesting  is  that  increasing  the  platform  inertia  did  not  seem  to  have  as  great  an 
impact  as  increasing  wheel  spin-axis  inertia. 

6.4  Simulation  Parameter  Effects 

The  simulation  parameters  also  affect  the  system  response  and  realism.  Although  one 
might  think  using  the  smallest  time  step  is  the  best  solution  (it  can  help  precision),  it 
is  well  known  that  decreasing  the  time  step  can  also  introduce  numerical  errors  [44] 
as  well  greatly  increase  the  run  time.  Figure  6.14  shows  the  results  of  different  test 
runs  with  different  fixed  integration  step  sizes. 

This  plot  shows  that  smaller  time  steps  are  more  accurate.  In  this  example, 
the  angular  momentum  magnitude  (which  is  required  to  be  constant  as  energy  is 
conserved)  is  compared  at  three  time  steps,  ts  =  0.1  sec,  ts  =  0.001  sec,  and  ts  = 
0.00001  sec  with  a  Runge-Kuttam  45  solution  method  (using  ode45  in  MATLAB) 
with  a  variable  step  size.  The  result  is  that  the  0.00001  sec  step  size  yields  an 
accuracy  in  angular  momentum  of  about  1.8  x  10^  —  6).  For  these  test  cases,  the 
angular  momentum  was  on  the  order  of  7470.0  kg  m2/  sec.  On  the  other  hand,  the 
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run  time  for  these  runs  was  greatly  increased  (for  10  seconds  of  simulation  time  run 
using  0.1  sec  took  less  than  5  minutes  whereas  a  0.2  second  run  of  the  0.00001  sec 
time  step  took  over  2  hours  to  run).  Thus,  a  simple  cost  benefit  analysis  shows  that 
there  is  a  balance  between  accuracy  and  run  time.  One  should  also  note  that  lower 
time  steps  do  not  always  yield  more  accuracy  as  round-off  errors  can  greatly  limit 
the  stability/performance  of  a  simulation  run.  So,  for  this  case,  a  time  step  of  0.1 
seconds  is  adequate  for  the  purpose  of  the  simulation. 

In  addition,  different  integration  methods  were  tested  using  variable  time  steps. 
As  far  as  using  ode45  (Runge-Kutta),  it  appears  to  only  affect  the  result  a  little  but 
wastes  much  more  state  storage  space  and  processing  time.  Thus,  for  this  case, 
Bogacki-Shampine  seems  to  be  the  best  method  to  use  in  terms  of  finding  a  middle 
ground  between  storage  space,  run  time,  and  numerical  accuracy /stability. 

6.5  Actuator  Test  Results 

The  last  primary  test  results  are  those  of  the  different  actuator  modeling  tests.  Fig¬ 
ure  6.15  is  a  power  actuated  versus  power  required  plot  combining  results  from  both 
the  case  of  modeling  the  actuators  as  identity  gain  and  as  uniform  scalar  gain.  The 
momentum  wheel  actuators  were  assumed  to  apply  torque  through  standard  DC  mo¬ 
tors  and  the  CMG  actuators  were  assumed  torque  the  gimbals  via  geared  DC  motors. 
Figure  6.16  shows  the  results  of  applying  three  actuators  aligned  along  the  B  axes 
instead  of  the  pyramid  configuration. 
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The  next  plot  shows  the  impact  of  power  tracking  when  using  the  different 
actuator  models.  Since  there  are  only  slight  changes  in  the  actuator  models,  as 
shown  in  the  preceeding  figures,  the  power  tracking  results  seem  to  be  very  close. 
But,  when  the  models  are  not  as  close  (e.g.  if  one  sets  the  scalar  gain  model  much 
lower  than  0.9),  then  the  result  is  a  not-so-perfect  power  output  track.  It  will  still 
track  the  desired  power,  just  not  as  quickly. 

Finally,  the  following  illustrates  the  use  of  three  actuators  instead  of  four  for 
the  simultaneous  attitude  control/energy  storage  system. 

One  can  see  from  this  result  that  the  system  performs  the  dual  functions 
marvelously  in  the  three  actuator  case.  However,  this  case  is  intended  to  verify 
the  flexibility  of  the  simulation  model  design  and  thus,  further  investigation  into 
using  three  actuators  versus  four  should  be  accomplished  in  terms  of  singularities, 
performance  degradation,  etc. 


6.6  Summary  of  Results 

This  chapter  presented  the  tests  results  of  those  tests  described  in  the  preceding 
chapter.  One  must  realize  that  a  large  number  of  simulation  tests  were  completed 
in  this  research  and  only  a  handful  of  these  tests  are  described  here.  Since  well  over 
100  different  simulations  were  run  during  this  research,  it  is  hard  to  include  plots 
of  all  the  effects.  The  only  important  effect  that  was  not  plotted  for  this  report 
relates  to  the  wheel  speeds.  It  should  be  noted  that  initial  conditions  seem  to  impact 
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the  simulation  results.  Namely,  the  initial  wheel  speeds  appear  to  impact  the  power 
tracking  profile  and  thus  the  simultaneous  solution.  It  was  observed  that  if  the  initial 
wheel  speeds  are  too  small  as  compared  to  the  peak  power  magnitude  and  duration, 
then  the  simulation  will  take  a  long  time  generating  the  required  power  and  torque 
when  the  power  profile  goes  through  a  peak  power  incident.  This  gives  some  slightly 
undesirable  results  in  the  simulation.  Even  though  the  required  power  is  tracked 
through  the  wheel  acceleration  and  not  the  wheel  speed,  the  wheel  speeds  do  have 
an  effect  on  the  total  output  power  magnitude  since  wheel  speeds  are  contained  in 
the  power  tracking  equation.  The  cases  where  the  norm  of  the  wheel  speed  initial 
conditions  were  small  was  that  there  seemed  to  be  a  slight  discontinuity  at  the  end  of 
the  peak  power  profile,  so  the  simulation  slowed  down.  This  problem  was  rectified  by 
increasing  the  initial  wheel  speed.  This  issue  seems  to  be  more  of  a  simulation  result 
and  should  not  be  a  problem  in  reality  as  the  wheels  will  work  in  a  normal  operating 
range  that  can  be  determined  from  the  anticipated  peak  power  requirement. 

In  order  to  recap,  the  results  here  were  broken  down  into  controller  comparison 
results,  system  parameter  results,  integration  method/time  step  results,  validation 
results,  and  actuator  modeling  tests.  From  these  tests,  one  can  gain  a  very  good 
understanding  of  the  effects  that  different  parameters  have  on  the  system  and  its 
control.  This  rewards  users  with  an  increased  ability  to  tailor  the  system  to  their 
own  needs. 
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Figure  6.1:  Regulation  Controller  I  Time  Response 
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Figure  6.2:  Regulation  Controller  II  Time  Response 
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Figure  6.3:  Regulation  Controller  Comparison  Vehicle  Time  Response 
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Figure  6.4:  Regulation  Controller  Comparison  Actuators  Time  Response 
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Figure  6.5:  Power  Tracking  Actual  Response  Versus  Profile 
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Figure  6.6:  Tracking  with  Controller  I  Time  Response 
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Figure  6.7:  Tracking  with  Controller  II  Time  Response 
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Figure  6.8:  Tracking  Controller  Comparison  Vehicle  Time  Response 
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Figure  6.9:  Tracking  Controller  Comparison  Actuators  Time  Response 
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Figure  6.10:  Control  Parameter  Vehicle. Time  Response  Test  Results 
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Figure  6.11:  Control  Parameter  Actuators  Time  Response  Test  Results 
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Figure  6.12:  System  Parameters  Vehicle  Time  Response 
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Figure  6.15:  Actuator  Power  Time  Response 
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Figure  6.16:  3-Actuator  Time  Response 
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CHAPTER  7 

CONCLUSIONS 


This  report  has  investigated  the  challenge  of  using  Variable  Speed  Control  Moment 
Gyroscopes  as  actuators  for  simultaneous  attitude  control  and  energy  storage  for 
spacecraft.  After  stating  the  problem,  the  relevant  literature  was  summarized.  This 
work  next  recapped  the  theoretical  development  of  an  algorithm  for  doing  simultane¬ 
ous  energy  storage  and  spacecraft  attitude  control  as  well  as  two  different  control  laws 
for  implementing  the  theory.  Then,  the  design  of  the  VSCMG  Workbench,  a  software 
tool  for  implementing  this  methodology  was  presented.  Fifth,  overall  system  and 
specific  controller  parameter  tests  for  analyzing  a  system  employing  this  system  on 
satellites  was  described.  Finally,  several  of  these  tests’  results  were  examined  in  order 
to  characterize  parameter  effects  for  an  example  simultaneous  attitude  control/energy 
storage  system. 

The  resulting  control  schemes  simultaneously  provide  attitude  control  torque 
and  energy  storage  torque  while  not  imparting  adverse  torques  on  the  satellite.  Al¬ 
though  both  controllers  respond  in  roughly  the  same  amount  of  time,  Controller  II 
seems  to  provide  a  much  more  desirable  result  as  the  vehicle  and  actuator  motion 


in  order  to  achieved  the  desired  attitude  are  much  smoother  than  controller  I.  As 


stated  in  chapter  6,  the  control  simulation  depends  greatly  on  the  actuator  initial 
conditions,  primarily  the  initial  wheel  speeds.  In  real  applications,  there  will  be  an 
operating  range  for  the  wheels  which  will  ensure  that  the  maximal  power  surge  to  be 
handled  by  the  satellite  during  eclipse  will  be  accommodated. 

The  test  results  indicate  that  some  parameters  affect  the  system  much  more 
than  others.  In  particular,  the  gimbal  inertias  do  not  seem  to  have  the  same  impact 
that  the  wheel  inertias  have,  especially  the  wheel  spin  axis  inertia.  Also,  it  seems  that 
the  most  relevant  controller  parameters  are  the  control  scalar  gain,  k,  and  the  matrix 
gain,  K.  The  mode  weights  do  not  seem  to  have  as  much  impact  as  the  gains. 
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CHAPTER  8 

RECOMMENDATIONS 


This  chapter  looks  at  several  different  areas  in  which  one  can  further  the  present 
research.  First,  one  should  look  into  characterizing  simultaneous  attitude  control  and 
energy  storage  singularity  avoidance  methods.  Related  to  this,  further  research  can 
examine  the  effects  of  wheel  speeds  on  the  resulting  vehicle  motion.  Specifically,  one 
can  examine  forcing  the  equalization  of  wheel  speeds  to  make  a  more  fault  tolerant 
system.  Equal  wheel  speed  solutions  would  allow  for  smoother  controller  operation 
once  a  wheel  fails  (e.g.  in  a  four  wheel  configuration). 

Second,  one  should  convert  the  user  defined  SIMULINK  S-functions  from 
MATLAB  *.m  files  to  C  code.  This  would  allow  faster  simulations  that  can  more 
accurately  be  used  with  Digital  Sequence  Processing  tools  for  real-time  implementa¬ 
tion  of  the  simulation. 

Third,  instituting  more  realistic  inertia  scalar  matrices  should  be  included  in 
the  simulation  model  and  perhaps  accounted  for  in  the  control  law.  This  includes 
both  the  overall  satellite  inertia  as  well  as  the  actuator  (wheel  and  gimbal)  inertias. 
This  would  allow  one  to  study  the  effects  of  flywheel  imbalance  as  well  as  allow  the 


control  designer  to  assess  the  validity  of  assuming  the  diagonality  of  the  wheel  inertia 
matrices. 

Fourth,  characterizing  different  disturbances  in  the  simulation  should  be  done 
as  well.  The  analysis  done  in  this  thesis  only  scratches  the  surface  as  far  as  disturbance 
rejection,  instead  primarily  focusing  on  the  problems  of  regulation  and  tracking.  One 
could  also  look  into  loop  shaping  and  other  performance  issues  such  that  the  control 
system  provides  the  best  system  performance. 

Fifth,  further  implementation  of  the  power  control  loop  which  includes  the 
power  bus,  its  voltage  regulation,  and  models  the  solar  array  power  conversion  method 
should  be  added  to  the  VSCMG  Workbench. 

Sixth,  one  needs  to  investigate  the  minimal  number  of  VSCMGs  required  for 
simultaneous  attitude  control  and  energy  storage.  This  would  relate  to  building  more 
fault  tolerant  combined  attitude  control/energy  storage  systems. 

Seventh,  much  more  fidelity  can  be  added  to  the  simulation  software  to  add 
realism.  More  realistic  torque  motor  transfer  functions  can  be  added  to  make  the 
system  respond  more  like  the  true  satellite  system.  In  this  vain,  better  models  of 
sensors  and  even  sensor  data  fusion  can  be  added  to  the  system.  The  simulation’s 
modularity  permits  this  model  fidelity  to  be  added  in  a  somewhat  straight-forward 
manner. 

Eighth,  many  different  nonlinear  control  laws  such  as  adaptive  control  laws, 
neural  networks,  and  robust  control  laws  can  be  can  be  added  to  the  system.  This 
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will  serve  two  purposes.  First,  it  give  one  more  insight  into  the  use  of  the  controllers 
in  order  to  see  which  control  laws  work  best  with  this  system.  Second,  this  will  help 
make  a  great  tool  for  inexperienced  control  system  designers  to  learn  the  effects  of 
these  control  laws  and  how  to  design  them  in  a  low-threat  environment. 

Ninth,  system  flexible  modes  can  be  modeled  in  the  simulation  as  well  as 
in  the  control.  The  system  developed  here  was  based  on  assuming  all  bodies  in  the 
system  are  rigid.  Flexible  analysis  will  permit  more  realistic  assessment  of  the  attitude 
control/energy  storage  system. 

Such  a  system  employing  the  derived  methodology  has  several  applications. 
Future  work  will  implement  this  control  methodology  on  a  more  realistic  space  vehicle 
simulation  model  to  be  used  for  controller  design,  hardware-in-the-loop  testing,  and 
vehicle  plant  modeling.  This  will  allow  a  low  cost  attitude  control  and  energy  storage 
system  design  prior  to  implementing  it  on  a  real  satellite.  Finally,  the  culmination 
of  this  work  will  be  a  ground  demonstration  of  flywheel  hardware  on  a  3  degree- 
of-freedom  spacecraft  simulator  at  AF.RL  to  validate  combined  attitude  control  and 
energy  storage  functionality. 
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APPENDIX  A 

DERIVATION  OF  THE 
ATTITUDE  DYNAMIC 
EQUATIONS 


In  order  to  derive  the  dynamic  equations  of  motion  presented  in  chapter  ch  3,  one 
should  first  analyze  the  motion  of  a  single  VSCMG  and  then  extend  the  results  to 
the  multiple  VSCMG  case.  This  is  done  by  deriving  the  angular  momentum  first  and 
then  finding  the  torque  acting  on  the  spacecraft  from  the  angular  momentum. 


A.l  Angular  Momentum 


Following  [42],  define  the  angular  momentum  for  one  VSCMG: 


h*ys/o  =  hP/o  +  hG/o  +  hw/o 


Using  the  parallel  axis  theorem: 


hP/O  —  hP/P *  +  hp*/o 
ho/o  =  ha/G *  _|_  hG*/o 

hw/°  =  Kw/w*  +  hw*/° 

Recall  that  W*  and  G*  coincide  by  assumption.  For  the  sake  of  completeness,  these 
values  have  been  left  in  the  equations.  This  distinction  is  dropped  later  in  this  deriva¬ 
tion. 

Generally,  as  mentioned  in  [42], 

-  A  V 

hX/°  = 

2  =  1 

So,  substitute  body  X’s  mass  and  treat  it  as  a  particle.  Then,  at  point  X *: 

hx"/ °  =  Yj  mi  {^ox  x  Nvx )  =  mx  (fox  x  Nvx ) 

2=1 

Nvx  =  (uJx/N  X  fox*) 

Assume  there’s  no  translational  motion  and  point  0  is  constant  in  A/". 

Now, 

hx  —  Tflx  (v ox*  X  CO x/N  X  7*  ox*) 


88 


Using  this  fact,  the  angular  momenta  of  the  particles  about  become 

hp*> o  =  mP  {fop*  x  Nvp *)  =  mP  (fop*  x  tiB/N  x  fop*) 
hG*/°  =  mG  (rog*  x  Nv° *)  =  mG  (rog*  x  ujB/N  x  fog*) 

hw*/o  =  mw  {fow,  x  Nvw *)  =  mw  {fow*  x  lub/n  x  fow *)  =  x  a;s/N  x  rog *) 

Since  P*,  G*,  and  W*  are  fixed  in  B,  one  can  see  that  the  angular  velocity  of  each  of 

the  3  points,  expressed  in  B,  is  cob/n- 

Next,  note  that 

hx*/0  =  mx  {fox*  x  ujx/n  x  rox*)  =  Ix*/°  •  £uy/n 

where  tiy/N  is  the  inertial  angular  velocity  of  point  -  mass  X*  and  /x*/°  is  the  inertia 
dyadic  for  a  point  mass  X*  with  respect  to  point  O,  as  used  by  [42].  This  is  defined 
as 

I X’/0  =  mx  (H(^x  •  rox*)  -  rox*fox*) 

Thus, 

^  "  ^Y/N  =  TnxifLf  ox*  '  p  ox* )  r  ox*  T  qX*  )  •  U)y/N  =  '^'X  (  f  ox*  '  T  ox*  )U  ‘  ^ Y/N  r  ox*  p  ox*  *  ^Y/n) 

F  °  '  & Y/N  =  'W'x{{y  ox*  '  T  ox*  )UY/N  T  ox*  fox*  ‘  W Y/N ))  =  ^x{f  ox*  ‘  pox*  )^Y/N  fox*  ‘  W Y/n)^ox *) 

F*/0  •  &Y/N  =  nT-xfox*  X  3y/N  X  -Vox*) 

Now, 

fix*/ O  =  |x*/o  .  tiY/N 
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Therefore, 


hpt/°  =  Ip*/°  •  ub/n 
Kg*/o  =  jg*/o  .  CoB/N 
h w*/o  =|w*/o  ,qb/n 

Furthermore,  the  Central  Angular  Momentum  of  a  body  Z  about  its  own  center  of 
mass  is  defined  as 

hz/z*  =  |z/z*  .  CJZ/N 


So, 


h?/ p*  =  Ip/p*  -cob/n 

ha/a*  =  p/ G*  •  Wg/AT 
fiw/w*  =  pv/w*  . 


Now,  calculate  the  system  angular  momentum  with  respect  to  point  O 

fisys/O  -  hP/O  _|_  ftG/O  4  pV/O  =  hP/P *  _|_  ftP*/0  +  /^G/G*  _|_  fr3*/0  4  ^W/W*  _|_  pV/O 

hsys/O  =  p/P*  .  $B/N  +  p*/G  .  t2B/N  +  p./O  .  UJ^i^r  +  pVVO  .  ^B/Ar  +  p/G*  .  OJqjn  +  |W/W*  . 

Re-grouping  terms  gives 

fisys/o  =  (|p/p*  +  p*/o  4  p*/o  4  pv*/o)  .  fiB/N  4  p/G*  .  fiG/N  4  |W/G*  .  Cjw/n 
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Next,  define 


JSC/O  =:  JP/P*  _|_  JP*/0  _|_  JG*/0  _|_  JW*/0 

Note  that  Isc/0  is  the  inertia  dyadic  of  all  inertia  that  remains  constant  in  B  over 
time,  it  includes  the  platform  inertia  with  respect  to  its  own  center  of  mass  plus  the 
point-mass  inertias  of  bodies  P,G,W  since  they  are  not  located  at  the  system  center 
of  mass. 

Rewriting  the  angular  momentum,  one  obtains 

h* wo  =  Is°/°  •  Cob/n  +  IG/G*  •  ug/n  +  Iw/G"  •  uw/n 

Next,  express  the  system  angular  momentum  about  point  O  in  terms  of  components 
in  B  as 

(h^°)B  =  [ISC/°]B  (Zb/n)b  +  [Ig/g‘]b  @o/n)b  +  [IW/W*]B  {*w/n)b 

Assume  gs,gt,9g  are  central  principal  axes  of  bodies  W  and  G  for  point  O.  Then 
the  inertia  matrices  of  body  G  about  G*  and  body  W  about  G*,  both  written  as 
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components  in  Q  are 


Igs  0  0 

[I<,'°1G=  0  IG,  0 

0  0  Io, 

Iws  0  0 

[Iw/g*]G=  0  IWt  0 

0  0  Iw, 

Now,  let  Lbg  represent  the  coordinate  transformation  from  Q  to  B  and  Lqb  = 
Ltbg  be  the  coordinate  transformation  from  B  to  Q.  Additionally,  per  the  angular 
velocity  addition  theorem 

&G/N  =  Wg/B  +  £>B/N 

and 
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Expressing  these  two  quantities  in  Q 


0 

(wg/b)g  = 

0 

7 

0 

ft 

(ww/g)g  = 

0 

(%/b)g  — 

0 

0 

7 

Next, 


(ojw/n)  b  =  Lbg  {p w/b)g  +  (3b/n)b 
Pg/n)  b  =  Lbg  {ug/b)g  +  (ub/n)b 


Thus, 


(ujw/n)  b  —  LBg  (<L>w/b)g  +  (<^b/n)b  =  LBg 


n 

0 


+  (&b/n)  j 


7 

0 


(3g/b)b  —  Lbg  Pg/b)g  +  (&b/n)b  =  T.bg 


+  ip  b  in') 


b 


7 


Now, 
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[I w'}]b  =  Lbg  [I °'g']gLtbg 


[Iw/g*]b  =  Lbg  [Iw/g1  glbq 

Here  note  two  (out  of  many)  different  ways  to  express  the  total  inertial  angular  mo¬ 
mentum  of  the  system  about  its  center  of  mass  written  in  terms  of  components  in  B, 
;  i)  Oh/Vidali’s  method  [32]  and  ii)  Schaub’s  method  [40]. 

Restate  the  angular  momentum  equation  just  derived 


(fcys/°)  =  [Isc/°]b  @b/n)b  +  [Ig/g1b  ($g/n)b  +  [Iw/w*]5  (uW/n)b 


This  can  also  be  written  as: 


\ 

0 

(h^°)B  =  [Isc/°]b  ($bin)b  +  Lbg  [Ig/g*]g  Ltbg 

Lbg 

0 

+  {ub/n)b 

7 

/ 

n 

\ 

+lbg  [i^]gltbg 

Lbg 

0 

+  (3b/n)  b 

7 

/ 
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Simplifying: 


0 

\ 

(to"/°)  =  [ISC/°]B  ( Ub/n)b  +  Lbg  po/o*]G 

0 

+  ^BG  (ub/n)  b 

7 

/ 

0 

\ 

+Lbg  [Iw/q1g 

0 

+  l^g  {3b/n)  b 

7 

/ 

Further  simplifying  yields 


B 


\ 

Igs 

0 

0 

0 

Mb  {3b/n) b  +  Lbg 

0 

IGt 

0 

0 

+  lTbg  (3b/n)b 

0 

0 

IGg . 

1 

7 

/ 

- 

" 

L 

\ 

Iws 

0 

0 

ft 

+Lbg 

0 

Iwt 

0 

0 

+  L'qq  (Cjb/n) b 

0 

0 

!wg 

7 

/ 

Now,  let 


9s  9t  9g 
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and 


(3b/n)g  =  Lqb  (3b/n)  b  = 


U!s 

gj  ( uB/n)b 

Ut 

= 

gj  {3b/n)b 

Ug 

_  $  (ub/n)b  _ 

Then  let  j  =  z\  +  z2  +  z%  where  z\  is  equal  to  [ISC/°]B  (cob/n) 

Deriving  z2 


B 


\ 

!gs 

0 

0 

0 

Z2  —  LBg 

0 

Io. 

0 

0 

+  Lbg  {ujb/n)  b 

0 

0 

!Gg  _ 

7 

/ 

9s  9t  9g 


Ig. 

0 

0 

0 

gj  {ub/n)  b 

\ 

0 

IGt 

0 

0 

+ 

9t  i^B in)  b 

0 

0 

Icg . 

7 

gTg  (Pb/n)  b 

/ 

^-Gs  9  s  ^Gt  9t  ^Gg  9g 


gj  {ub/n)b 

sT  (< 1b/n)b 

7  +  gj  (pB/N)i 

=  Igs gsgj  {ub/n)b  +  Ig tgtgf  (# b/n)b  +  Ig eg9gJ  {^b/n)b  +ihgg9 


<B 
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Next,  Z3  becomes: 


z3  =  L 


BG 


\ 

\ 

Iws 

0 

0 

Cl 

0 

Iwt 

0 

0 

+  Ltbg  (wB/n)  b 

0 

0 

!wg 

\v 

/ 

(gS  9t  9g^j 


Iws 

0 

\ 

0 

Cl 

gj  {3b/n)b 

\ 

0 

Iwt 

0 

0 

+ 

gj  (3b/n)b 

1° 

0 

Twg ) 

7 

9g  (Pb/n)b 

/ 

Z3  = 


Iws  9s  Iwt  9t  Iwg  9 9 


^  +  9J  (w b/n)b 
gj  (<2b/n)b 
i  +  gj  ($b/n)b 
~  (<2b/n)  b  +  Iwt 9t9f  {^b/n)b+  Iwg 9g  9g  (Pb/n) 


B 


+OIwsfl's  +  T^Wg  9g 

From  this  one  can  write  the  angular  momentum  equation  for  one  VSCMG  in 
a  different  way  [40] 


(hsys/°'j B  =  [ISC/°]B  (ojb/n) b  +  (Igs+Iws)  9s9j  (<^b/n)b  +  (Ig(  +  Iwt)  9tgJ  b/n ) B 

+  (lGg  +  Iwg)  9ggJ  {ub/n)b 
+7  (lGg  +  Iwg)  9g  +  ^Iws  9s 
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A.  2  Torque 


Now,  find  the  total  external  torque  of  the  system  about  point  0  with  one  VSCMG: 


fsys/o  _  NjL  ^sys/o^  =  Bd_  ^SyS/o^  _j_  (Qb/n')b  X  h^0 
Next,  break  up  this  value  into  2  terms.  The  first  term  is  given  by: 


B  d 


7,  (*"'°)  =  Bi  (Ms  (“b/»)b)  +  Bi  Mis  (aa/»)B  +  “ft  Mis  (‘V/s) 


\  ) 


B 


and  the  second  term  is  calculated  as 


(&B/n)  b  X  fiSyS//°  =  OJb/N  [Isc]B  (^B/n)  b  +  Ub/nLbG 


(. 


Igs  o  0 
0  IGt  0 

0  0  In 


/,  „ 

Iq.  0  0 


+3b/nLbg 


0  IGt  0 


V 


0  0  Ic 


lbg^b/n  +  Wb/nLbG 


\ u  u  jg9j  yrj 

(  \  ( 

I\vs  0  0  Cl 

0 


7 


0  Iwt  0 


0  0/, 


Wa 


Iws  0  0 


\V 


+3B/nLbg 


\ 


0  Iwt  0 


0  0  Iwa 


rjp  ^ 

Lbg^b/n 
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Returning  to  the  first  term,  one  can  substitute  in  the  angular  velocity  terms: 


0 

\ 

\ 

Bi  '  [Isc,o]B  (Jb/k)b  +  BA 

LBG[IW'}GLTBG 

Lbg 

0 

+  (^b/n)b 

7 

) 

I 

( 

( 

n 

\ 

\ 

1  B  d 
^  dt 

Lbg  [Iw/G‘](j  Ltbg 

Lbg 

0 

+  (ub/n)b 

7 

) 

/ 

The  above  equation  reduces  to: 


I 

O 

_ i 

\ 

Bi  (S*"'°)  =  Bi  [iso/°1b  (°b,k)b  +  Bi 

Lbg  [Ig'g*]g 

0 

V 

1 

_ 1 

/ 

0 

\ 

\B  d 
"r  dt 

Lbg  [Iw/g‘]g 

0 

V 

7 

/ 

+  (Leg  {^'q-]0LIg 


+  B3(lbg{I°'°-}gLtbo 


Next,  take  the  derivatives  (note:  a  single  dot  over  a  vector  epresents  its  time  derivative 


taken  in  frame  B.  A  single  dot  over  a  scalar  represents,  of  course,  its  scalar  time 


(wB/jv)J 
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derivative) : 
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Now,  substitute  the  inertia  scalar  matrices  into  the  equation  for  IG/Q*  and  Iw/°* 


Bd_ 

dt 


(  hsys/°  \  =  [Isc/o] 


B 


+ 


Lbg 


+ 


Lbg 


+ 


■'BG 


/ 


+ 


Lbg 


( 


+ 


JBG 


( 


+ 


L 


BG 


h3 

0 

0 

Igs 

0 


( '^B/N ) 


B 


+ 


0  0 

iGt  o 
o  i0$ 
0  0 

IGt  0 


0  0  IGq 


Ig.  0 


0  I( 


0 

Gt  0 


0  0  IGo 


0 

0 

7 


\ 


f 

“ 

\ 

Igs 

0 

0 

0 

jsc/o 

B  B/N )  B  + 

Lbg 

0 

0 

0 

0 

0 

lag 

7 

y 

\ 

Ig,  0 

0 

0 

+ 


jy 


-'BG 


\ 


L'bg  (3b/n) 


b 


( 


+ 


0  lot  0 

0  0  Ia 

Igs  0 

o  i* Gt 


Lbg 


lbg  {3b/n) 


b 


+ 


Lbg 


0 

7 

0 

0 


0  0  7, 


Ga 


\ 


LTbg  (3b/n) 


b 


Js  0  0 


0  Jt  0 


Iws 

0 


0  0 
Iwt  0 


0  0  IwQ 


Iws 

0 


0  0 
hvt  0 


0  0  Iwa 


I\Vs 

0 


0  0 
iwt  o 


0  0  Iwa 


) 

[ 

0 

0  J 

\ 

Cl 

Iws 

0 

0 

0 

+ 

Lbg 

0 

I\vt 

0 

7 

0 

0 

1 _ 

\ 

ci 

Iws 

0 

0 

0 

+ 

Lbg 

0  . 

Iw, 

0 

7 

) 

0 

0 

IWg 

\ 

7 

7 ws 

0 

rT  f  \ 

JBG  [UB/N) B 

+ 

Lbg 

0 

7 wt 

) 

0 

0 

\ 


L^g  {ub/n} 


b 


n 

o 

7 


\ 


jy 


\ 


L^  (vb/n) 


B 


lW„ 


\ 


Lbg  {ub/n) 


b 


101 


Next,  we  use  fact  that  Isc/°,  I°/G* ,  and  Iw/G*  are  constant  with  respect  to  B,  Q,  and 
G,  respectively.  Thus,  the  listed  inertia  derivative  terms  are  zero,  so: 


'  Ig.  0  0  0 

ft  (^sys/°)  =  [Isc/°]b  b  +  Leg  0  IGt  0  0 

i  0  0  IGg  7 


f  Igs  0  0  0 

+  Leg  0  IGt  0  0 


Ig .  0  0 


Lbg  0  IGt  0  lbg  (3b/n)j 


0  0  IGg  [7 

Igs  0  0 


0  0  IGa 


Igs  0  0  Iws  0  0  Q, 

+  Lbg  0  IGt  0  LTBG  (uB/n)  b  +  LBG  o  IWt  0  0 


_  o  0 

Iws  0  0  1  [  tl 


0  0  Iw  7 


Iw .  0  0 


+  Leg  0  IWl  0  0  +  LBG  q  IWt  0  LTBG  {ojB/N) 


0  0  Iw  7 


0  0  h 


Iws  0  0 

+  Lbg  o  IWt  o 


Js  0  0 


L'eg  ft b/n )  B  +  LBg  0  Jt  0  LTBG  (tiB/N^J 


0  0  Iwa 


0  0  Jn 


Substituting  the  fact  that  LBG  =  LBG 7  and  =  —AfL'BG,  where: 


0  -(7)  0 
^=700 


(A.1) 


0  0  0 


102 


and  also  noting  that: 


B  d_ 
dt 


IGs  0  0 

7  0  IGt  0 

0  0  IGg 

then  we  are  left  with: 


^sys/oj  _  [Isc/o]b  ^  + 


0 

0 

-7 

0 

0 

0 

0 

— 

7 

0 

0 

0 

— 

0 

7 

0 

0 

0 

IGg  7  _ 

0 

/ 


-'BG 


( 


IGs  0  0 
0  IGt  0 


0  0  IGg 


J 

1 

o 

\ 

0 

i _ 

/ 

(A-2) 


‘ 

' 

\ 

/ 

' 

' 

\ 

Ig3 

0 

0 

Igs 

0 

0 

+ 

Lbg{  7) 

0 

0 

LBg  (£b/n)  b 

+ 

Lbg 

0 

ht 

0 

(“7 )L'bG  B/N )  B 

0 

0 

IGg_ 

) 

V 

0 

0 

JGg_ 

/ 

+ 


Lbg 


lGs 


0  0 


0  I, 


Gt 


\ 


L BG  B/N ) 


B 


+ 


Abg(t) 


lws 


0  0 


0  L 


wt 


0 


\ 


I 

0 

0 

IGg\ 

)  \ 

\ 

[ 

0  0 

Iwg\ 

v\ 

\ 

I\vs 

0 

0 

0 

/ 

ws 

0 

0 

+ 

Lbg 

0 

I\vt 

0 

0 

+ 

Lbg  (7) 

0 

Iwt  0 

L'bo  (ub/n)  b 

0 

0 

IWg 

7 

) 

0 

0 

Iwg 

) 

\ 

I\vs 

0 

0 

Iws 

0  0 

+ 

Lbg 

0 

Lvt 

0 

(“7 )L'bq  (3b/n)  b 

+ 

L 

BG 

0  / 

wt  0 

L'bg  ( 

*b/n)b 

0 

0 

Lw, 

) 

I 

0 

O 

1 _ 

/ 
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Combining  term  1  and  term  2  yields: 


■fsys/O  =  NJt  (hsys/°)  =  Bj-t  +  (ub/n)b  X  ^sys/0  =  [ISG/°]b 


B 


' 

■ 

\ 

( 

r 

\ 

Ig. 

0 

0 

0 

Ig.  0 

0 

+ 

Lbg 

0 

lGt 

0 

0 

+ 

Lbg{  7) 

0  IGt 

0 

LTbg  (3b/n)  b 

0 

0 

7 

) 

0  0 

lG°. 

J 

\ 

Igs 

0 

0 

1 

1 

Ig. 

0  0 

+ 

Lbg 

0 

iGt 

0 

( 

-7)^bg  {“b/n)b 

+ 

Lbg 

0 

lot  0 

Lbg  ( 

0 

0 

Io9_ 

) 

I 

0 

0  IGg_ 

\ 

/ 

Iws 

0 

0 

n 

+ 

Lbg(  7) 

0 

Iwt 

0 

0 

+ 

0 

0 

IWg 

7 

/ 

\ 

Iw. 

0 

0 

+ 

Lbg{  7) 

0 

Iwt 

0 

L'bg  b/n )  B 

0 

0 

Iwg 

\ 


B 


+ 


Lbg 


Iwa  0  0 

0  I\Vt  0 

0  0  Iw„ 


\ 


L BG  (mb/n'J 


B 


' 

‘ 

‘ 

> 

Iw. 

0 

0 

n 

Lbg 

0  1 

wt 

0 

0 

0 

0 

IWg 

7 

J 

\ 

7 

-1 

- 

\ 

Iws 

0 

0 

+ 

Lbg 

0 

Iwt 

0 

( 

“7 )Lbg  (to b/n)  b 

/ 

\ 

0 

0 

<0 

_JI - 

J 

Ig. 

0 

0 

0 

+  ub/nLbg 

0 

Ig* 

0 

0 

0 

0 

- , 

Ol 

7 

Js 

0 

0 

O 

O 

4" 

q 

+lob/n-Lbg 

0 

Jt 

0 

LbqGjb/n  +  to  b/n  Lbg 

O 

4 

O 

0 

0 

0 

J9 

1 

O 

O 

7 

+  &B/N  [I  SC  1b  (&b/n) 


B 


Note  that  the  above  equation  refers  to  one  VSCMG.  We  can  extend  this  equation 


104 


to  the  multiple  n-VSCMG  case  by  summing  all  the  contributions  from  the  gimbal 
structures  and  wheels.  If  we  define: 


[Itot/o]s  =  [lsc/o]s  + 


/  n 


(££**  ([iGi/cnGi + [iWi/onGi)  hQiBt 


and  substitute  1°/°*,  and  Iw/G*  back  into  the  equations,  this  can  be  written  more 
compactly  as: 


(7®>/°)b  =  [/T°T/°]fl  (cSb,n)  +  ( 

“  - 

\ 

+ 

EWr«i0, 

i=l 

0 

7 i 

/ 

’ B 


\ 


+  LffiB  (Pb/n)  j 


/ 


+ 


B 


Ei®o1([i,vof]0l  +  r',of]0i) 

i—  1 


V 


o 

o 

7 i 


\ 


J  J 


+  -  L*°<  ([I0,,°;]g.  +  r,,0f]ft).(il)io.*)  («»/»), 


0 

\ 

0 

\ 

a 

\ 

+ 

ELBGi  [IGi/G*]Gi 

i=l 

0 

+ 

z— 1 

0 

+ 

z=l 

0 

7 i 

J 

7< 

J 

0 

/ 

\ 

+  (<3b/N)b 

n 

£W  lw^]Gi 

0 

0 

/ 
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A. 3  Simplifying  the  Torque  Equation 


Next,  let 


7i 

1 

\ 

\ 

Oi 

Dx 

72 

A 

Ew^r^'ic, 

2=1 

0 

+  ■f'GiB  {ub/n)  b 

7i 

) 

/ 

7n 

7i 

0 

\ 

I>2 

72 

-  (wb/»)s 

EiB0,([i0''0r]ft  +  [iw,,0;k) 
2=1 

0 

7* 

/ 

7i 


D* 


72 


-  (  E  LBGiiii)  [IGi/anGi  LGiB  -  Lbg, 

<i= 1 


In 


7i 

0 

\ 

72 

A 

e  iBGi  [i°',anGi 
2=1 

0 

7n 

7* 

/ 

7i 

0 

\ 

72 

A 

±LBGir^}Gi 

2=1 

0 

7n 

7* 

/ 
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fil 

fi2 

Q>n 


\ 

n 

££BCirron0, 

i=l 

0 

0 

/ 

fix 

fi2 

Cln 


cf 

\ 

tLBGi  [iWi/°nGi 

0 

i=  1 

1 

o 

1 _ 

/ 

This  yields  the  system  dynamic  motion  equations  (presented  chapter  3 


(fsys/o) ^  =  [Itot,o]b  (iiB/K)  +  (ob,n)  (3b/n)b 

i—  — i  r  *i  r  n 


7i 

7i 

fix 

+  ( Di  +  -D2  +  D%) 

72 

+  (#1  +  B2) 

72 

+  £ 

fi2 

+  F 

7n 

7n 

fin 
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APPENDIX  B 

SIMPLIFICATION  OF  THE 
ATTITUDE  DYNAMICS  - 
N OTATION AL  DERIVATIONS 


This  appendix  presents  the  derived  simplification  of  the  dynamic  equation  of  motion 
presented  in  chapter  3.  In  order  to  derive  this  result,  one  must  first  review  the  relevant 
matrices  liberally  used  in  this  derivation. 


G gm  - 


9gi 


9gn 


9gi  9gi  9gi 


Ggd  =  blockdiag 


i 


LL 


9gi 


9gn  9gn  9gn 
■  9gn 


!Ggk  ~ 


^Gg  1  •••  I Ggn 


IGgd=dmg{lGgk} 


These  hold  for  the  spin  and  transverse  axes  of  the  wheels  and  gimbals  as  well. 
Since 

Jcjg  =  leg  +  hvg 

one  can  define  the  following 

JGgd  =  hgd  +  Iwgd 

u>d  =  blockdiag  |  J^b/at)  (&b/n)  (ub/n)  ’  ’  *  (^b/w)  ]  | 
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Now,  derive  explicit  expressions  for  the  B\,  B2,E,  F,  D\,  D2,  and  Dz  matrices: 


The  derivation  is  identical  for  B2,  E,  and  F  except  one  uses  Iwgd,  Iwad,  Gs,  and  Iw3d 
plus  Gs  pre-multiplied  by  ui  to  find  the  respective  matrices.  Furthermore,  D2  is  found 
similarly  to  jB4 
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The  matrices  derived  are  summarized  as 

Bi  =  G9IGgd  = 

fel®lGgl)  (5fl2®lGg2)  ••• 

fen®Icgn) 

B2  =  Gglwgd  = 

fei®lwgl)  fe2  ®  iwg2 )  •• 

■  fen®  I  Wgn) 

II 

£ 

CO 

O 

II 

(9 si  <s>  iw.i )  (gs2 

®  Iw,2  )  •  •  • 

fen  ®  W.J 

F  —  (ub/n'J BGsIwad  —  B 

(9s  1  ®  Iw51 ) 

fe2  ®  IwS2 ) 

D2  =  (fiB/N^j  ^GgJgd  =  (cOB/n') 

ST 

<0 

C-| 

OP 

( 9 92  ®  ^g2  ) 

•  fen  ®  IW,„  ) 

fen  ®  Jgn) 


111 


112 


APPENDIX  C 

CONTROL  LAW  DERIVATION 

DETAILS 


C.l  Controller  I 

In  this  section,  support  details  in  the  derivation  of  the  control  law  for  controller  I 
(derived  in  Chapter  3)  is  presented.  Note  that  this  control  law  is  very  similar  to  that 
used  in  references  [37],  [32],  and  [39].  It’s  main  assumptions  are  that  one  is  interested 
in  regulation  of  the  vehicle  about  some  fixed  inertial  reference  and  can  thus  neglect 
orbital  motion.  Since  most  of  the  control  law  is  derived  in  Chapter  3,  only  details 
not  included  earlier  are  described  here. 

First,  note  that 

Cl  (C.l) 

and 


Vi  =  ke[e  i 


(C.2) 


In  addition,  if  one  defines 


where 


and 


results  in 


Then, 


Q(P)  = 


~Pl  —&2  Pz 
Po  Pz  @2 
Pz  Po  ~Pl 
—@2  Pi  Po 


Po 

Pi 

P2 

Pz 


QWr)  = 


Pro 

Pn 
Pr  2 
Pr3 


~Pn 

1 

to 

-Pro 

Pro 

1 

$*> 

CO 

Pn 

Pn 

Pro 

-Pn 

~Pr2 

Pn 

pro 

QT  (Pr)  P  =  ~QT  (P)  Pr 


(C.3) 


(C.4) 


(C.5) 


(C.6) 


(C.7) 
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This  is  true  since  the  left  hand  side  of  the  above  equation  is  equal  to: 


—fin  fid  +  fir o  fil  +  firzfi 2  “  fir2fi 3 
~  fir  2  fid  ~  fir  2  fil  +  fir  o02  +  finfiz 
—  fir  2  fid  +  fir2fil  ~  fin  fit  +  firofiz 


(C.8) 


Now,  one  computes  the  right  hand  side: 

1  fir,  r  1 

—fil  fid  fiz  —fi>2  —filfiro  +  fid  fin  +  fizfir2  —  fi2fir3 

fin 

— /?2  fiz  fid  fil  ~  ~  —fi2firo  ~  fiz  fin  +  fid  fir  2  +  fil  fin 

-  fin 

—fiz  fi2  —fil  fid  —fizfiro  +  fi2fin  —  fil  fir  2  +  fidfirs 

fin 

Slightly  rearranging  this  result  yields: 

—fid  fin  +  filfiro  —  fizfir2  +  fhfin 

fi2fir0  +  fiz  fin  —  fid  fir  2  ~  fil  fir  z  (C.10) 

fizfiro  ~  fiifin  +  fil  fir  2  ~  fidfirs 

Now,  comparing  the  left  hand  side  to  the  right  hand  side  verifies  that: 

QT  (fir)  fi  =  ~QT  (fi)  fir  (C.ll) 

These  relationships  are  useful  in  the  subsequent  derivation.  One  should  also  note 
that  since  efei  is  a  scalar,  then  its  transpose  is  equal  to  efei.  This  relationship  will 
also  be  needed  here. 
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C.1.1  Lyapunov  Equation  Analysis 


Realize  that 


Simplifying  this  yields 


Vi  =  ke\e  i  +  ke^ii 


Vi  =  2 ke\Te\  =  2k  (f3T  -  ftr)  (/3  -  ft) 


Multiplying  this  term  out  becomes 

Vi  =  2 A;  (fp  -ffP-  PTPr  +  ftTft) 

Next,  note  that: 

ftrft  =  ||ft||2  =  1 

and 

PTP  =  \\pf  =  l 

Then,  taking  the  derivatives  of  these  terms  yields: 

PTP  =  PTP  =  PjPr  =  PtT  Pr  =  0 


(C.12) 

(C.13) 

(C-14) 

(C.15) 

(C.16) 

(C.17) 


Now,  the  V\  equation  is: 

Vi  =  2k  (-ftr/3  -  £rft)  =  2k  -  £rft)  (C.18) 

As  in  [32],  the  Euler  parameter  formulation  yields  the  kinematic  differential 
equation 

P  =  \q(P)u  (C.19) 
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and  similarly 


A  =  \q{Pt)ut  (C.20) 

Now,  plug  this  into  the  V\  equation  to  get 

Vi  =  2k  ^-^ujjQT(pr)/3  -  ^nrQT(/3)A^  (C.21) 

plugging  the  earlier  relationship  that  QT  ((3r)  /3  =  —QT  (/?)  / 3r ,  yields 

Vi  =  2k  (^ujQT(p)pr  -  (C.22) 

Finally,  this  becomes: 

Vi  =  -  (ooT  -  ujJ)  {kQT(m)  (C.23) 

This  completes  the  details  that  supplement  the  text  in  chapter  3. 

C.2  Controller  II 

In  several  applications  requiring  three-axis  attitude  control,  devices  such  as  Reaction 
Wheel  Assemblies  and  Control  Moment  Gyroscopes  must  not  only  provide  the  neces¬ 
sary  torques  but  also  precision  is  required.  For  this  reason,  a  slightly  more  accurate 
control  law  is  investigated  here  to  yield  different  results  for  the  simultaneous  problem. 
The  theory  derived  is  similar  to  the  one  presented  earlier  [8]. 

First,  as  shown  in  [45],  5Cj  can  be  simplified.  Use  the  strap-down  equation  to 
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show  that  Lbr  —  as  follows 


Lrn  —  —&Lbn 

(C.24) 

Lrn  =  —&tLrn 

(C.25) 

LBr  =  LBnLrn  +  LbNLrn 

(C.26) 

=  LBN(—tirLRN)T  H - (<^Lbn)L^n 

(C.27) 

=  LbnL^n(—( Cjt)t  H — (do)LBNLjiN 

(C.28) 

=  LBR(—&r)T  H - ( U))LBR 

(C.29) 

=  LBR&r  H — &Lbr 

(C.30) 

LbR^t  =  LbrOJt^t  —  O^Lbr^t 

(C.31) 

=  —GjLbr^t 

(C.32) 

—  — <jj(oj  —  Scj) 

(C.33) 

=  w&j 

(C.34) 

Second,  one  can  simplify  Soj 

Sdj  =  do  —  LsROJr  —  dj5uj 

(C.35) 

=  d)  —  LbrUIt  —  a )(cd  —  LbrUt) 

(C.36) 

—  do  —  LbrOJt  +  dbLsROJr 

(C.37) 

(C.38) 

Note  that  this  relationships  will  become  important  as  the  derivation  unfolds. 
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C.2.1  Approach 


Prior  to  solving  for  the  control  law,  one  can  summarize  the  approach. 

1.  At  time  t,  one  is  given  given:  j3r,  u>r,  ojt,  j3 ,  u>,  7,  and  7. 

2.  Solve  for  rr  to  yield  stable  regulation/tracking  (this  solution  is  derived  in  the 
remainder  of  this  Appendix). 

3.  Using  Tr,  one  finds  Cl  and  7  to  provide  stable  attitude  control  and  power  track¬ 
ing. 

4.  Assuming  less  than  perfect  actuators,  one  actuates  7  and  Cl  through  Tome  and 
rwiric,  the  motor  command  torques. 

5.  Use  the  motor  output  torques  tgvtio  and  Twmo  to  find  the  actuated  7  and  Cl. 

6.  Integrate  to  get  the  new  7  and  Q. 

7.  Repeat  step  1  at  t  —  t  +  At. 

C.2.2  Stable  Attitude  Controller 

Begin  with  the  Lyapunov  function  candidate  (where  k  >  0): 

V  =  kVq  +  Vu  (C.39) 

where  Vq  =  S(5\  +  5/3%  +  +  (5f3 0  —  l)2  and  Vw  — 
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One  can  see  that  this  Lyapunov  candidate  is  positive  definite  and  radially 


unbounded.  Thus,  it  is  a  valid  Lyapunov  function  candidate.  Next,  find  V. 


V  =  kVq  +  Vu  (C.40) 

Vq  =  2{50\50\  +  50 2502  +  50^50^  +  2(500  —  1 )  ( <5  /5o ) )  (C.41) 

50o 

r  i  5k 

Vq  =  2  5(3 o  5(3 1  5/3 2  -^o  (C.42) 

L  J 

H  =  2(50T50  -  50o)  (C.43) 

Also,  one  finds  VL  as  follows: 

K,  =  \5ujtIt5uj  +  \5u)tIt5lo  +  \5(jjt'It5ijJ  (C.44) 

Note  that  \5u)tIt5uj  =  \5ljtIt5u,  so  our  equation  is  now 

Vw  =  5ujtIt5lo  +  ^5uitIt5oj  (C.45) 

Thus, 

V  =  2k(50T50  —  5/3o )  +  5cotIt56j  +  15ujtIt5lo  (C.46) 

Page  418  of  reference  [8]  gives 

50 1 

25P0  =  —5ut  s/32  (C.47) 

50 3 
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This  implies  that 


-60  o  =  \6ujt 


60 1 
602 
60s 


(C.48) 


Also, 


And  this  implies  that 


60  =  \Q(60)8u 


60t  =  \5ujtQt  (80)  (C.50) 

Next,  substitute  these  relationships  into  the  V  equation,  yielding 


(C.49) 


+  6cutIt6uj  +  | 6lotIt6uj  (C.51) 


f 

\ 

60 1 

V  =  2  k 

\6u)tQt  (60)60  +  \6ojt 

60 2 

60 3 

/ 

Then, 


V  =  k6ujTQT  (60)60  +  k6u i1 


60 1 
60 2 
60s 


+  6(jJTIt6Co  +  h6uT  It6<jo  (C.52) 


And, 


L 

1 

r  -| 

\ 

60 1 

h 

3 

<-o 

II 

kQT  (60)60  +  k 

602 

T  It66j  T  \It6w 

60s 

/ 

(C.53) 
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Now,  define  a  positive  definite  gain  matrix,  K  so  that: 


—K5u)  =  kQT(8(3)5/3  +  k  +  ItSoj  + 


(C.54) 


From  this,  one  can  see  that 


V  =  —8cutK8lu  <  0 


(C.55) 


This  implies  that  the  system  is  Lyapunov  Stable  since  V  is  less  than  or  equal  to  zero. 
In  order  to  draw  a  conclusion  about  asymptotic  stability,  one  must  employ  LaSalle’s 
theorem  by  examining  the  largest  invariant  set  contained  within  the  set  of  all  states 
where  V  is  zero.  Clearly,  since  K  is  positive  definite,  then  for  V  to  be  zero,  Soo  must 
be  zero.  That  is  to  say 


-KSoj  =  kQT(5(3)5/3  +  k  +  ItSuj  +  \It8oj  =  0 


(C.56) 


This  further  implies  that 


kQ^  (<5 /3)5 0  +  k  +  ItSui  —  0 


(C.57) 


Additionally,  within  the  largest  invariant  set  where  V  is  zero  (which  in  turn 


implies  that  5co  is  zero)  necessitates  that  5uj  is  zero  for  all  time.  Obviously,  for  the 
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case  where  5u  is  zero  for  all  time,  then  most  assuredly  5u>  must  also  be  zero  for  all 
time  in  order  for  6cu  to  be  contained  in  the  largest  invariant  set. 

One  must  now  identify  where  V,  6u>,  and  5u  equal  zero.  Here,  the  stability 
constraint  addressed  in  (C.57)  reduces  to 

50x 

kQT(60)60  +  k  S02  =0  (C.58) 

60s 

This  can  be  simplified  even  further  when  one  examines  QT  {5 0)50: 

i  t  f  -i 

—5f3i  — 5(32  5/33  5/30 

S/3 o  -5/33  5f32  6/3i 

'  (C.59) 

5/33  5/?o  —5/3 1  5f3  2 

—S(32  5f30  5/3  3 

This  results  in 

— 5/?i  5/3q  +  S/3qS/3i  +  5/335/32  —  Sf328f3^ 

-60250 o  -  +  60o502  +  50x50z  =  0  (C.60) 

—50^500  +  6025(3]  —  60x602  +  50o50s 
Now,  one  is  left  with  a  case  where 

60x 

k  60 2  =  0  (C.61) 

60  3 
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This  implies  that  5  Pi  =  <5/52  =  5P$  =  0.  Since  5/3  must  have  a  norm  of  1,  then 
given  that  three  of  the  4x1  column  matrices  are  zero,  5Po  must  equal  1.  So,  the 
largest  invariant  set  where  V  is  zero  contains 

1 
0 

(C.62) 

0 
0 

From  this,  one  can  see  that  5/3  goes  to  the  identity  transformation  over  time, 
which  implies  that  the  orientation  difference  from  the  current  angular  position  and  the 
desired  angular  position  is  equal  to  identity.  Thus,  the  body  is  at  its  desired  angular 
position.  To  summarize  the  stability  result,  when  equation  (C.54)  is  satisfied,  then 
the  largest  invariant  set  where  the  positive  definite,  radially  unbounded  function  V’s 
derivative  is  zero,  the  only  states  contained  in  this  set  are  at  the  origin  (i.e.  the  points 
where  the  attitude  position,  velocity,  and  acceleration  errors  are  zero).  This  implies 
global  asymptotic  stability  of  the  closed  loop  system. 

Now  that  stability  has  been  proven,  one  only  now  needs  to  find  controls  that 
satisfy  (C.54).  Recalling  that 

55)  —  CO  —  LrrUIj-  +  CcLbrUj-  (C.63) 
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One  needs 


\ 

K5uj  =  — 

&Qt(<W£  +  A; 

602 

+  It(6j  —  LbroJt  +  6jLbrujt)  +  | It6u> 

60z 

/ 

(C.64) 


This  can  be  rearranged  to 


K(u  —  Lbr^t)  =  —It^>  +  ItLbr^t  ~  It&LbrWt  —  ^It6lo  —  k 


60 1 

602 

Sfc 


(C.65) 


Now  note  that 


It  ui  =  tb  —  (D/j’ca  —  By  —  Dsy  —  EVl  —  FVl 


(C.66) 


This  yields 


Ku>  —  KLbr^t  ~  ~ Td  T  CjIt^J  +  By  +  Dsy  +  Etl  +  FQ,  +  ItLbru>t  (C.67) 

60 1 

I 

ItuLbrWt  —  \It{u  —  Lbr^t)  ~  k 


60 2 
60  3 


(C.68) 


Further  rearranging,  this  becomes 


Koj  —  KLsRWr  T  td  —  CjItw  —  ItI'BR^t  T  It^jLbr^t  —  k 


60 1 
602 
60 3 


=  —\It{u  ~  Lbrut) 


+By  +  Dsy  +  EQ,  +  F(0!.69) 
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Now,  let 


This  is  similar  to  the  control  law  derived  earlier.  Therefore,  the  required  torque  is 
summarized  as 

50 1 

tr  =  Kto  —  KLbrwt  —  CoItu  —  ItLbr^t  +  IrSoLbr^t  ~  k  (C.73) 

S/3 3 

It  should  be  pointed  out  that  once  one  obtains  the  attitude  error  8/3,  then  one 
can  get  the  transformation  Lbr  (from  1Z  to  B)  using 

(50l  +  5/31  -  5(31  ~  8 Pi)  2(80,802  +  80,50,)  2 (iftift  -  m) 

=  2(8 0,5 02  ~  50,50,)  (50l  -  50\  +  801  -  801)  2(502803  +  50,500) 

2(50,80 3  +  W&)  2(50,80 ,  -  WA)  (5/302  -  <5/32  -  «5/32  +  8 ft) 

(C.74) 
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So,  the  question  remains,  how  does  one  compute  5P 1,  6(3 2,  and  5/33  at  any 
instant  in  time?  From  page  443  of  reference  [8],  one  finds  that 

5  Pi  (3o  Pz  p2  ~Pi  Plr 

5  P2  ~Pz  Po  Pi  —P2  Plr 

(C.75) 

<%  P2  ~Pl  Po  Pz  Pzr 
&Pi  Pi  P‘2  Pz  Po  Por 

Therefore,  one  is  able  to  find  the  attitude  error  to  compute  the  required  torque  from 
equation  (C.73)  for  stable  tracking  of  the  vehicle  reference  attitude. 
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APPENDIX  D 

ACTUATOR  CONFIGURATION 


Figure  D.l:  Generic  VSCMG  Layout 


Figure  D.2:  VSCMG  Pyramid  Configuration 
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APPENDIX  E 

EXAMPLE  COMBINED 
ATTITUDE /POWER 
TRACKING  PROBLEM 


Most  of  the  orbital  example  is  derived  from  Tsiotras,  Shen,  and  Hall  [5,  4].  This 
orbital  example  of  a  satellite,  Iridium  (spacecraft  number  25578),  provides  a  rigorous 
case  study  in  the  simultaneous  attitude/power  tracking  control  system  under  worst 
case  polar  orbit  conditions.  In  this  scenario,  the  flywheel  energy  storage  function 
augments  a  satellite’s  solar  array  system.  During  sunlight,  the  solar  array  provides 
energy  to  the  many  spacecraft  subsystems.  Any  excess  power  is  used  to  charge  the 
flywheels  within  the  VSCMG  suite.  Then,  during  eclipse,  the  satellite  uses  the  array 
to  power  its  many  subsystems.  Since  there  is  no  back-up  for  the  VSCMGs  during 
eclipse,  the  system  must  be  prepared  to  meet  a  peak  power  requirement. 

Note  that  a  few  RWA  control  laws  for  this  problem  were  developed  in  ([45]). 


Figures  (E.l)  and  (E.2)  reflect  the  satellite  with  solar  arrays  deployed,  the 


tracking  scenario,  the  associated  coordinate  reference  frames. 

I 
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Figure  E.2:  Scenario  Layout 
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APPENDIX  F 

VSCMG  WORKBENCH: 
SIMULATION  DESIGN  DETAILS 


Iscvec 

Iscl 


n  Inertia  Overview 


34 


Figure  F.3:  Internal  System  State  Derivative  Computation  Logic 


Figure  F.4:  Gimbal  Rate  Torque  Amplification  Constraint 


Figure  F.5:  Actual  Plant  Output  States  Oscilloscope 
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FINAL  NOTE 


“Man’s  flight  through  life  is  sustained  by  the  power  of  his  knowledge.”  (Austin  ‘Dusty’ 
Miller,  the  quote  on  the  Eagle  and  Fledgling  Statue  at  the  US  Air  Force  Academy). 


